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Abstract 


In this thesis we consider a system of two linear partial differential equations of parabolic type, 
which represent the basic behavior of temperature and vapor density with the effect of evapo- 
ration and are coupled in a particular way. In the first part, we prove the existence and unique- 
ness of the solution of this system, first in a domain of one spatial dimension and then in a 
domain in R° delimited by two horizontal planes. To do this, we construct a particular variant 
of Fourier series. In the second part we consider a system of analogous equations in a spherical 
domain. It is the system of linear parabolic equations representing the behavior of the tem- 
perature and the density of vapor in the case of a spherical droplet with evaporation. Using a 
variant of Fourier series relative to a spherical domain, we prove the existence and uniqueness 


of the solution of this system in the case of spherical symmetry. 


Key-words: System of linear parabolic equations, Fourier series, discontinuous coefficient, 


case of spherical symmetry, evaporation. 
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Résumé 


Dans la présente thèse nous considérons un système de deux équations aux dérivées partielles 
de type parabolique linéaires représentant le comportement essentiel de la température et de 
la densité de vapeur avec l’effet de l’évaporation couplées d’une manière particulière. Dans la 
première partie nous démontrons l’existence et l’unicité de la solution de ce système d’abord 
dans un domaine d’une dimension spatiale et puis dans un domaine dans R* délimité par 
deux plans horizontaux. Pour ce faire nous construisons une variante particulière de la série 
de Fourier. Dans la deuxième partie nous considérons un système d'équations analogues dans 
un domaine sphérique. Il s’agit du système d'équations paraboliques linéaires représentant le 
comportement de la température et de la densité de vapeur dans le cas d’une gouttelette sphé- 
rique avec l’évaporation. En utilisant une variante de la série de Fourier relative au domaine 
sphérique, nous démontrons l’existence et l’unicité de la solution de ce système dans le cas de 


symétrie sphérique. 


Mots-Clés : Système d'équations paraboliques linéaires, série de Fourier, coefficients disconti- 


nus, cas de symétrie sphérique, évaporation. 
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Introduction 


d f E thème de la présente thèse est l’étude d'équations linéaires aux dérivées par- 
tielles de type parabolique, couplées d’une manière particulière. Le choix de 
ce thème est motivé par l'intérêt pour la description mathématique du phénomène 
d’évaporation de l’eau à partir de la surface d’eau, et, plus en général, des interactions 
entre l'air et l’eau. En effet, l’évaporation de l’eau est un des facteurs principaux de l’in- 
teraction entre l’air et l’eau dans notre environnement, mais sa description mathéma- 
tique rencontre beaucoup d'obstacles techniques. On peut estimer que ces difficultés 
résultent fondamentalement du fait que l’évaporation de H:0 se produit sur la surface 
de l’eau, c’est-à-dire sur une variété de dimension 2, et crée par la chaleur latente une 
source négative de la chaleur concentrée sur une variété de dimension 2, tandis que 
l'on est intéressé ses conséquences dans le domaine de dimension 3 occupé par l'air 
ou par l’eau. 

Au début de notre recherche, nous cherchions une modélisation des interactions 
générales entre l’air et l’eau. Mais les circonstances particulières autour du phéno- 
mène d’évaporation nous ont conduits à l’étude d’un modèle simplifié de la densité 
de la vapeur d’eau et de la température avec l’évaporation. C’est un modèle simplifié, 
car les difficultés mathématiques intrinsèques qui résultent de la singularité du phé- 
nomène physique nous ont obligés à considérer deux équations linéaires décrivant, 


dans leur approximation linéaire, la température et la densité de la vapeur d’eau avec 
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leur diffusion. 

Malgré la description approximative par un modèle simplifié, le problème formulé 
dans le cadre de l’Analyse mathématique présente - nous croyons — des aspects inté- 
ressants qui méritent notre attention. Le couplage de deux équations est formulé par 
une source de la chaleur concentrée sur une variété de dimension inférieure à celle du 
domaine et par une condition sur la frontière du domaine. Pour surmonter les diffi- 
cultés dues à cette situation particulière nous avons introduit et analysé une série de 
Fourier particulière, ce qui constitue le point crucial de notre travail du point de vue 
technique. 

Pour avoir une idée un peu concrète sur le système d'équations que nous allons 
analyser dans la présente thèse, citons ici le cas le plus simple de notre système d’équa- 
tions. En effet, pour une fonction T représentant la température dans l’eau (—b < x3 < 
0) et dans l'air (0 < x3 < a) et une fonction II représentant la densité de la vapeur dans 


l'air (0 < x3 < a), on considère le système d'équations 
CoodtT = 0x3 (KO x3 T) + WO(X3) pour {>0,-b<x3<a, (1) 


OH=Yo0,,I1 pour 1z0,0< x3 <a, (2) 
avec les conditions de couplage 


II(#,€1) — II(E, 0) 


€1 


y) = Yi (3) 


[| (4) 


=7T0+@T| 


X3=0 X3=0? 


et d’autres conditions usuelles sur les frontières et des conditions initiales (Ô(x3) dans 


(1) est la delta de Dirac sur R 3 x3). La relation (3) est une approximation de la condition 


: OI (5) 
ven Ôx3 x3=0” 
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tandis que (4) est une approximation de la condition 
[| ,, 0 = Tos(T)|x, 0 (6) 


Tys(T) étant la densité de la vapeur saturée à la température T (voir (1.11) dans le cha- 
pitre 1). Dans (2)-(4) Yo, Y1, «1 sont des constantes strictement positives, tandis que Cw 
et x dans (1) sont des constantes strictement positives séparément dans —b < x3 <0et 
dans 0 < x3 < a. Même si (1) et (2) sont des équations du type parabolique linéaires, 
les conditions de couplage (3)-(4) et le fait que c,, et x ne sont globalement pas des 
constants nous obligent à trouver une méthode particulière pour résoudre ce pro- 


blème. 


Contenu de la thèse 


La présente thèse est composée de quatre chapitres. Dans le premier chapitre, nous 
allons illustrer la motivation physique et mathématique de l'interaction entre l'air et 
l’eau ainsi que les méthodes principales utilisées. Plus précisément, nous allons rap- 
peler le système d'équations mécanique et thermodynamique qui décrivent les inter- 
actions de l'air et l’eau, y compris le phénomène de l’évaporation de l’eau. Ensuite 
nous allons illustrer l’idée des méthodes principales que nous allons utiliser dans la 
présente thèse. 

Dans le deuxième chapitre, nous allons considérer le système de deux équations 
paraboliques linéaires, dont la formulation est motivée par une modélisation approxi- 
mative de la température et de la densité de vapeur dans l’air et dans l’eau avec l'effet 
de l’évaporation qui se produit sur la surface de l’eau. Dans l’équation modélisant la 
variation de la température, on trouve les coefficients de la diffusion de la chaleur et la 
chaleur spécifique différents dans l’eau et dans l’air ainsi que la source (négative) de 
la chaleur concentrée sur l’intersurface entre l’eau et l'air. Pour traiter cette équation 


nous introduisons une version particulière de la série de Fourier. En outre, en utilisant 
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cette série de Fourier, nous allons démontrer l'existence et l’unicité de la solution sta- 
tionnaire de ce système d'équations dans un domaine d’une dimension spatiale. 

Dans le troisième chapitre, nous allons considérer premièrement le système d’équa- 
tions linéaires de la température et de la densité de vapeur dans le domaine réduit à 
un domaine de dimension 1, c’est-à-dire sur l'intervalle ] — b, a[ CR et nous allons dé- 
montrer l'existence et l’unicité de ce système. Puis nous allons démontrer l'existence et 
l’unicité de la solution du système complet dans le domaine Q = {(x1, x2, x3) € R° | — b < 
x3 < a} sous une hypothèse restrictive (voir (3.1) dans le chapitre 3), en utilisant la va- 
riante de la série de Fourier introduite dans le deuxième chapitre et aussi des espaces 
de Sobolev définis par la série de Fourier pour x3 et la transformée de Fourier pour x 
et x». 

Dans le quatrième chapitre nous avons un problème analogue dans le cas sphé- 
rique. Il s’agit du système d'équations paraboliques linéaires représentant la tempéra- 
ture et la densité de vapeur dans un domaine sphérique {x € R°||x| < b} avec le do- 
maine occupé par l’eau {x € R°||x| < a}, 0 < a < b, comme dans le cas d’une gouttelette 
d’eau sphérique avec l'effet de l’évaporation. Nous allons introduire une série de Fou- 
rier relative au domaine sphérique avec la symétrie sphérique, suivant l’idée de la série 
de Fourier introduite dans le chapitre 2, mais avec une élaboration assez consistante 
pour l’adapter au domaine sphérique. En utilisant ces outils, nous allons démontrer 


l'existence et l’unicité de la solution dans le cas de la symétrie sphérique. 
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1.1 Motivation physique et mathématique : Interactions 


air-eau 


6 f ES interactions entre l’air et l’eau interviennent dans beaucoup de phénomènes 
C2) quiintéressent notre vie. Nous pouvons citer des exemples de phénomènes 
violents comme cyclone tropical ou “el Niño”, pour lesquels l'interaction entre l’atmo- 
sphère et la mer est un facteur essentiel, mais aussi le climat local déterminé par la 
présence d’un lac ou la proximité de la mer, jusqu'aux phénomènes à une échelle as- 
sez petite comme la surface du thé versé dans une tasse ou celle d’une gouttelette d’eau 


de la pluie. On trouve une littérature abondante pour expliquer, ou tenter d'expliquer, 
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ces phénomènes comme par exemple [8], [21] pour de grandes échelles et [12], [36] 
pour de petites échelles. 

La surface de l’eau exposée à l’air est, comme nous le savons bien, capable de pro- 
duire l’évaporation selon les conditions physiques (pour les fondements physiques du 
phénomène d’évaporation, voir par exemple [20]). Or, l’air ainsi que l’eau sont des 
fluides suscetibles d’être mus et ceci, comme nous l’observons communément, pro- 
voque la variation de la position de la surface d’eau. Donc pour décrire d’une manière 
générale les interactions entre l’eau et l'air, nous avons besoin des principes de la ther- 
modynamique qui déterminent les conditions de l’évaporation de H,0 ainsi que des 
équations de la mécanique des fluides (pour les principes de la thermodynamique, voir 
par exemple [24]; pour les équations de la mécanique des fluides, voir par exemple 
(251). 

Dans chaque partie — la partie de l’air et la partie de l’eau - considérée séparément, 
le mouvement du fluide et la distribution de la température sont décrits par un sys- 
tème d'équations. Dans la partie de l’air on peut considérer par exemple le système 


d'équations 


Ôd#(b+7)+V:-((p+7)v)=0, (1.1) 
Ôfn + V-(xv) = AA, (1.2) 
ul CE 
@+mGrv+v-Vu)=nAv+( + )V(V-v) RvILe Fe JT}-(p+xmgez, (13) 
a h 
@+mMaôT +v-VT)+R{(E + T)T}v.v= (1.4) 
Ha  Hnh 


3 Ôv; ÔVE 2 Ô 
= K1AT + + "6 ;$Vv) —v;+0(V:v)? + En. 
K] 12 Gx ôx; 3 jk Dr &(V:v) 1 


Où 0, x, v, T désignent respectivement la densité de l’air sec, la densité de la vapeur, la 
vitesse, la température, et 1, 1, &, C1, K1 sont respectivement le coefficient de diffusion 
de la vapeur dans l'air, le coefficient de viscosité d'écoulement, le coefficient de visco- 
sité volumique, la chaleur spécifique de l'air et le coefficient de diffusion de la chaleur 


dans l’air, tandis que -gez = —-£(0,0, 1)? est la force gravitationnelle par unité de masse 
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et FE est l’éventuelle source de la chaleur; ici la pression est donnée par 


p=R( +7, (1.5) 


Ha  Hh 


La et Un étant respectivement la masse molaire de l’air sec et celle de H,0. 
D'autre part, dans la partie occupée par l’eau liquide, le mouvement de l’eau li- 


quide sera régi par les équations de Navier-Stokes 
dsv+(v-V)v+Vp=vAv+ges, (1.6) 


V:v=0, (1.7) 
et la fonction T représentant la température devra satisfaire à l'équation 


3, Ov; Oôvx, Ô 
Cd T+v-VT=KkAT+v Ÿ (+ 


— )— v;+ E. 1.8 
‘ei Ôxk 6x) 8x d 8) 


Dans (1.6) et (1.8) v, ©2, k2, E2 sont respectivement le coefficient de viscosité de l’eau 
liquide, la chaleur spécifique de l’eau liquide, le coefficient de diffusion de la chaleur 
dans l’eau liquide et l’éventuelle source de la chaleur. 

Le système d'équations (1.1)-(1.4) est assez compliqué et, même si le système de 
Navier-Stokes (1.6)-(1.7) est bien connu et bien étudié, le système d'équations (1.6)- 
(1.8) lui aussi n’est pas simple. Mais dans la modélisation mathématique de l’interac- 
tion entre l’air et l’eau ce qui rend particulièrement difficile le problème est la présence 
de l’inter-surface, qui est en principe inconnue (donc les deux sous-domaines occupés 
par l’air et par l’eau liquide sont eux aussi en principe inconnus). 

Sur l’inter-surface, que nous notons $, nous devons imposer avant tout la contrainte 
géométrique : la vitesse du déplacement de S dans la direction de la normale ñn à S doit 
être égale à 


-n, (1.9) 


où v( et v@) sont la vitesse de l’air et celle de l’eau liquide. 
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Deuxièmement la continuité du tenseur de contrainte, complétée par la tension super- 


ficielle, exige, dans le cas où l’intersurface S peut être définie par une fonction h(x1, x), 
S — {X3 = h(x, X2)}, 


l'égalité 


- 2 
= pOn;= [niv +0 ;vl)) +(C- DV: ya UCI +0;v®)]n;+ (1.10) 


i i 


j=1 


Vh 


Ve )n,, 
TENTE UE 


+ 


où p) est la pression de l'air donnée par (1.5) et p®) est la pression de l’eau liquide 
(qui est notée p dans (1.6), tandis que y est le coefficient de la tension superficielle 
et le terme yV: RS est l'effet de la tension superficielle (pour les détails de la 
tension superficielle, voir le Chap. VII, $ 60 de [25]). 

En ce qui concerne l’évaporation de l’eau, qui se produit sur la surface d’eau, son 
aspect fondamental peut être expliqué dans le cadre de la physique statistique (voir 
[20], [3], [36], etc). Mais, pour que l’on puisse analyser l'effet de l’évaporation dans 
l'air et dans l’eau, on a besoin d'introduire les relations macroscopiques concernant 
l’évaporation. Pour cela, nous rappelons d’abord la densité de la vapeur saturée, que 
nous notons 7,,(T) ; elle est en effet essentiellement fonction de la température T. Sa 


valeur établie par les physiciens est 


_ l 7,63(T-273,15) 
Tys(1)=——EÆ0-10 TAHX , Ep =6,107 (mbar), (1.11) 
RIT 
où 
R 
R = — 
Un 


(nous avons pris cette formule de [27], où se trouve aussi l’explication de la formule 


(1.11)). Ceci signifie que dans l’air à la température T la densité de la vapeur d’eau 
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ne peut pas dépasser la valeur x,.,(T) ; dans le cas où à cause de la diminution de la 
température 7,,(1) devient inférieur à la densité réelle de la vapeur d’eau, la partie 
excédante de la vapeur doit se condenser et devenir liquide (ou solide) en formant par 
exemple un brouillard. D'autre part, le principe de l’équilibre local sur lequel se base 
la description macroscopique implique que dans le voisinage immédiat de la surface 
de l’eau la densité de la vapeur doit être égale à celle de la vapeur saturée. 

L'autre quantité essentielle pour le phénomène d’évaporation est la chaleur latente 
de la transition de phase de l’eau de l’état liquide à l’état gazeux (et vice-versa). Sa 
valeur est 


Lir(T) = (8244 — 2,72 T)10° U1kg). (1.12) 


Quand se produit l’'évaporation de l’eau, ceci prive l’énergie thermique L;,(T) par unité 
de masse de l’eau évaporée. La chaleur latente L;,(T) elle aussi dépend de la tempéra- 
ture T, mais sa dépendance de T est relativement petite, à différence de la densité de 
la vapeur saturée 7,,(T), dont la dépendance de la température T est très significative. 

Le problème du mouvement de deux fluides avec une intersurface dans le cadre de 
la mécanique, c’est-à-dire avec les conditions (1.9) et (1.10) mais sans les conditions 
dues à l’évaporation, a été étudié par certains chercheurs depuis le travail de Tani [40]. 
Le problème de la surface libre d’un fluide visqueux incompressible (c’est-à-dire, pro- 
blème réduit au seul côté de l’eau) a été étudié, depuis les travaux de Solonnikov [38], 
[39], par plusieurs auteurs [16], [15], [41]. 

Sur l’évaporation de l’eau à partir de la surface de l’eau on trouve de nombreux 
travaux de caractère physique mathématiques ([12], [31], [3], [29], [19], [6], [36], [43},, 
etc.). Mais, comme le problème de l’intersurface entre l’eau et l’air, même si on le 
considère seulement dans le cadre de la mécanique, est assez difficile, il nous semble 
que pour le moment l’étude dans le cadre de l’Analyse mathématique du système 
d'équations décrivant le mouvement des deux fluides avec l’intersurface et l’effet de 
l’évaporation est difficilement réalisable. 


Compte tenu de ces circonstances se pose la question : comment on peut décrire 
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le phénomène de l’évaporation à partir de la surface de l’eau et ses effets — la diffusion 
de la vapeur d’eau dans l'air et la diffusion dans l’air et dans l’eau de l’effet thermique 
de l’évaporation — sous l'hypothèse de l’absence du mouvement macroscopique de 
l'air et de l’eau. Même dans cette simplification des conditions, il nous semble que 
l'analyse satisfaisante du système d'équations du modèle qui décrit ce phénomène 
n'est pas facilement réalisable. Comme nous l’avons évoqué dans l’Introduction, à 
notre avis, ces difficultés sont dues au fait que les deux équations sont coulpées par les 
conditions sur l’intersurface comme (5), (6). Notre intérêt principal réside donc dans 
la compréhension de la possibilité de la description de ce problème par un système 
d'équations aux dérivées partielles et de ses propriétés à partir de l'existence et l’uni- 


cité de la solution. 


1.2 Méthodes principales 


1.2.1 Équations paraboliques linéaires 


Les équations aux dérivées partielles du type parabolique et les méthodes de leur 
résolution sont bien connues; on peut les trouver facilement dans les manuels des 
Équations aux dérivées partielles ([28], [11], etc..). Pour l'intérêt de la problématique 
de la présente thèse, nous rappelons la méthode de la résolution de l'équation de la 
chaleur dans le domaine d’une dimension spatiale à l’aide de la série de Fourier; pour 
simplifier, choisissons le domaine 


1=]0,7{ 


et considérons l'équation de la chaleur 


2 


Ô 
—U(É,X) =K 


: oz UE, + f(x) dans [0,oo1 x] (1.13) 
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avec les conditions aux limites 
u(t,0) = u(t,7) =0 (1.14) 


et la condition initiale 


u(0, x) = uo(X) dans 1, (1.15) 


où K,Q est une constante strictement positive et w,(x) et f(f, x) sont des fonctions don- 
nées. 


On rappelle que les fonctions 


2 
2e k=1,2,:.:, (1.16) 
VT 


constituent une base orthonormale de L°(0, x) et un système orthogonal complet de 
2 .. _ 
H*(0, x) et sont des vecteurs propres de l'opérateur — L avec la condition de Dirichlet 
homogène (c’est-à-dire annulation de la fonction aux points x = 0 et x = x), vecteurs 
propres correspondants aux valeurs propres À4£ = k?, c’est-à-dire 
d? v2 v2 
—sinkx=k—=sinkx,  kK=1,2,.... (1.17) 


dx? 7 VT 
Donc, si on pose 


uok= Î LE sin xd ñ@= | Fu, 0 sin kxdx k=1,2,--.-, (1.18) 


on peut décomposer l’équation (1.13) avec les conditions (1.14)-(1.15) en une famille 


de problèmes de Cauchy 
d 2 
PAU = —Kok"ux (D + (D, ug(0) = uox  K=1,2,.... (1.19) 


La solution u(f,x) du problème (1.13)-(1.15) sera alors construite par les solutions 
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uy.(t) des problèmes de Cauchy (1.19) dans la forme 


O0 v2 
u(É,X) = uy(f) — sin Kx. 
2407 


Nous voulons appliquer l’idée de cette méthode à l'équation 
v(OÔU(E, x) = d,(k(HOutr, ») +y(D6(%) dans R:x]-— Li, Lol, 
avec les conditions aux limites 
u(t,—Li) = u(t, Li) =0 


et la condition initiale 


U(0, X) = up(X). 


Ici les coefficients v et x sont tels que 
V(x)=v1 si —Li<xX<0, V(xX)=v2 si 0<x<L), 


K(X)=K1 Si —L<x<0, Kk(x)=Kk2 si0<x<L», 


V1, V2, K1, K2 étant des constantes strictement positives. 


Nous considérons aussi l'équation 
v(r)ôçu = 4, (r”k(r)ô;u) + yô(r — a) dans R; x ]0,bf, 


avec 
v(r)=v si0<r<a, V(r)=v sia<r<b, 


K(T)=Kk] Ssi0<r<a, K(r)=Kk2 sia<r<b, 


V1, V2, K1, K2 étant des constantes strictement positives. 


(1.20) 


(1.21) 


(1.22) 


(1.23) 


(1.24) 


(1.25) 


(1.26) 


(1.27) 


(1.28) 
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L'équation (1.26) sera envisagée avec la condition initiale 
u(0,r) = uo(r) pour re]0,bl, (1.29) 
et la seule condition aux limites 
u(t,b) =0 pour {2>0; (1.30) 


au point r = 0 on ne donne pas de condition aux limites, ce qui est conforme au fait 
que le domaine ]0, b[ n’est autre que le domaine {x € R°|0< ,/x° + x? + x? < b} avec la 
symétrie sphérique des fonctions considérées. 

Les équations (1.21) et (1.26) ne sont pas équivalentes à l’équation (1.13). Mais si 
on peut trouver une série de Fourier alternative, alors on pourrait construire la solution 
des équations (1.21) et (1.26) dans une forme analogue à (1.20). 

Les autres informations sur l'équation de la chaleur et celle de la diffusion de masse 


qui peuvent être utiles peuvent être trouvées dans [2], [30], [42], etc. 


1.2.2 Espaces de Hilbert appropriés 


Soit Q un domaine borné de R” muni de la frontière régulière 0Q. Si (a; ;(x));,j=1,.,n 


est une matrice carrée réelle d'ordre ñn définie pour tout x € Q telle que 


aij(x)= aji(x) VG,j) el... ,n}x{l,...,n} (1.31) 
D ajéié;zxlél  VéeR”,Kk>0, (1.32) 
i,j=1 


alors, comme il est bien connu, on peut définir le produit scalaire 


à s Fr Lu Ôv(x) 


dx = (u, v}f 1.33 
Oj,j=1 ôx; Ôx; (u, V)fn(o) (1.33) 
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et avec ceci l’espace de Hilbert 
H'(Q) = {u : mesurable, Q — R|{u, u) <o, u = 0 sur 00}, (1.34) 


dont la norme est 


ul fnço, = Qu, u»)?. 
Les vecteurs propres eK de l'opérateur elleptique 


n 


Ô Ô 
(Au) (x = - Y es a (Gi uto), (1.35) 
J 


i,j=1 


c'est-à-dire les fonctions ex, k = 1,2,--., vérifiant la relation 
Aex = Àrer (1.36) 


avec les nombres À£ > 0 (les valeurs propres de l'opérateur A), constituent une base 
orthonormale de L?(Q) et un système orthogonal complet de H'(Q). La décomposition 
orthogonale d’une fonction f € L?(Q) et l'application de la relation (1.36) nous permet 


de résoudre immédiatement le problème 


- 2 Le (ao | u(x)) = f(D, (1.37) 
0x 
j=1 Xj 


& 
Il 
© 


sur 00. (1.38) 


Pour les détails de cette méthode, voir [28], Chapitre IV, Section 1. 
Une fois construite cette base orthonormale {ex}. de L’(Q), en utilisant cette base 
orthonormale et la décomposition orthogonale, on peut résoudre sans difficulté non 


seulement l'équation elliptique (1.37), mais aussi l'équation parabolique 


n 


Ô Ô Ô 
au À 8x (a na x)) = f(x), (1.39) 
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avec la condition aux limites 
u(t,x) = 0 pour (f,x)eR; x 0Q (1.40) 


et la condition initiale 


u(0, x) = up(x) pour xeQ (1.41) 


(voir par exemple [28], Chapitre VI, Section 2, Point 2). 
Pour appliquer cette idée à nos équations (1.21) et (1.26), il est essentiel de choisir 


les espaces de Hilbert appropriés. Nous définissons les espaces de Hilbert 
L (1 ={u:1—R, mesurables | [rcoiucorax < œo!} (1.42) 
I 
muni du produit scalaire et de la norme 


Ga vam [ vœux Iulz2çn = 4/ (u U)r2(n 


et 


H\(D={ueC(TR) IÉCTORE < oo, u(—Li) = u(Li) =0} (1.43) 
I 


muni du produit scalaire et de la norme 


(ut =|IrxœurDdx, lalan=/mGüa. 
OT], 1(D 1(D 


Nous définissons également 


b 
L?(0, b) = {@:(0,b) —R, mesurables | Î v(r)r{lu(r) dr < co}, (1.44) 
0 


b 
A} (0,b)={peL2(0,b) l kr Eu dr < co, u(b) = 0}, (1.45) 
0 
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munis respectivement du produit scalaire 


b b 
GA, v)2n = [ v(r)r£u(r)v(r)dr, (U,V)f0,b = [ k(r)r£u/(r)v'(r)dr, (1.46) 


et munis de la norme correspondante à ces produits scalaires. 

Comme on le verra, le couple d’espaces de Hilbert L2(1) et H}(1) définis dans (1.42) 
et (1.43) comme le couple d'espaces de Hilbert L2(0, b) et H}(0, b) constitue la base de 
la structure dans laquelle nous allons résoudre nos équations. 


Nous utiliserons aussi les espaces de Hilbert 


b 
L°(a, b) = {@ : (a, b) —R, mesurables | Î r’lu(r)l?dr <oo}, (1.47) 
a 


b 
H"(a, b) = (pe Tab | AE u(nÉdr < 00, u(a) = u(b) = 0} (1.48) 


munis respectivement du produit scalaire 
b b 
(U,V) tab sl ru(r)v(ndr, (uv) | ru(r)v'(r)dr (1.49) 
a a 


et munis de la norme correspondante à ces produits scalaires. 


1.2.3 Variantes de série de Fourier 


La méthode principale sur laquelle nous nous appuyons pour traiter les équations 
de la chaleur (1.21) et (1.26) est l’utilisation d’une version particulière de série de Fou- 
rier. 

Les équations (1.21) et (1.26) ne sont pas identiques à l'équation (1.13). Mais elles 
ont des aspects similaires à ceux de l’équation (1.13), de sorte que nous pouvons uti- 
liser une méthode similaire à la résolution du problème (1.13)-(1.15) par la série de 


Fourier (1.16). 
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Plus précisément, on définit la famille de fonctions eg(x), k= 1,2,---, par 


yk(x) 


er (x) = , 
lykllz2cn 


Agv 


yk(x) = 


FU sin (BU + V/ 2x) 


. pe , 
L sin | LL +4/ ui x) si -L1<x<0, 
V'kivi Ki # 


(1.50) 


(1.51) 
si 0O<x<l). 


Ici Ag, BU»), YU) sont choisis de telle sorte que y£(x) s’annule au point x = L, et que 


la continuité de la fonction y£(x) ainsi que celle de (x) & yk(x) soient garanties au 


point x = 0 (pour la définition précise de À+, B (&), YU»), voir (2.16), (2.17), (2.20) dans 


le chapitre 2). 


Comme il sera illustré dans la Proposition 2.3.1, la famille de fonction {e,}®, sera 


une base orthonormale de l’espace de Hilbert L£(1) et aussi un système orthogonal 


complet de H}(1). Il est important que les fonctions eg(x), k = 1,2,--, jouissent de 


5 


bonnes propriétés similaires aux fonctions = sin &x, qui forme la série de Fourier clas- 


sique sur l'intervalle ]0,7[. On verra en effet que la Proposition 2.3.1 garantit les pro- 


priétés suivantes : 


: 2 


d d 
ii) 7 EU) et) = Ag v(x)er(x) 


iii) K(X) Lex(x) est continue sur [—L;,L)]; 


iv) sup 


—Li<x<L) 


lex(x)l < K1, 


VxEe] — Li,0[01]0, Lol 
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1 1 K K 
Ki = 2max| | }mex[ ve ve) 
VViK1 VV2Kk2 LVL 
T 2 
v) Mok-1)<A;<Motk+1),  Mo= | —— >) | 
L — Vi 2 


De manière analogue, pour le cas sphérique on peut construire la base orthonor- 
oe) 2 


male {ek}®, de L2 (0, b), qui est également un système orthogonal complet de A} (0, b). 


Plus précisément, on définit 


er(r) = ui , (1.52) 
Ye (1720, 
où 
sin (B(Ax)+ ML (an) 
YA —— 5 pour 0<r<a, 
yk(r) = n (1.53) 


sin ( KE? (b- n)) 
EE our a<r<b 
Agv2 Tr P mn: 


Ici Àk est une suite réelle strictement croissante pour tout k e N, B(Ax) et y(A+) 
sont choisis de telle sorte que la continuité de la fonction yz(r) ainsi que celle de 
K(r) _ y&(r) soient garanties au point r = a (pour la définition précise de À+£ voir (4.51), 
et pour celle de B(A+), y(+), voir (4.32), (4.33) (avec 1 = #), dans le chapitre 4). 

Alors, outre que {e4}® , est une base orthonormale de L2 (0, b) et un système ortho- 
gonal complet de H!(0, b), on verra, dans la Proposition 4.4.1, que la famille de fonc- 


tions {e4}®. jouit des propriétés suivantes : 
k=1 


: 2 
i er|l< = À4; 
) Vertou = À& 


ii) ek(r)et r2x(r)-Lex(r) sont continues ; 
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a d 
iii) ——er(r)lr=0 = 0,  ex(b)=0; 
dr 
: d 2 d 2 
iv) Tr K(r)——er(r) = Akcwo(r)r"eg(r) Vre]0,a[0]a,bl; 
dr dr 


v) ilexiste une constante C1 < c telle que 


lek(a)l < Ci VkeN\{0}, 


vi) eg(r) possède k — 1 zéros sur ]0, bl; 


VKIK2 5 


SR M pen 


Ces propriétés seront utilisées pour démontrer l'existence d’une solution des sys- 


tèmes d'équations que nous considérons. 


1.2.4 Espaces de Sobolev anisotropes 


Dans le cas du problème dans le domaine de dimension 3, plus précisément dans 
le domaine Q = R°x]— b, af, nous allons construire la solution dans un espace de type 
espace de Sobolev avec des ordres différents dans les directions de x1, x2 et dans la 
direction de x3. 

Plus précisément, nous allons utiliser l’espace de Hilbert H°"(Q), s > r, défini par 
la série de Fourier pour x; et la transformée de Fourier pour x1 et x, caractérisé par la 
norme 


1/2 


(ul sr ço = | > [ Q+kéP) "(++)" IA, dé) (1.54) 
k=1 
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ici A(É, k) est la fonction de ë € R? et de k € N\{0} définie par 
U(, K) = F(xi,x2) (ER U)12(-b,a) 


avec la transformée de Fourier (%,,x)(+) par rapport à (x1,x2) et le produit scalaire 
(er, U)12(-p,a dans l’espace de Hilbert L?(-b, a), ex étant un des éléments de la base 
orthonormale {e;}®, de l’espace L?(-b, a) définie dans (1.50). Même si les fonctions 
eg ne sont pas des fonctions trigonométriques, comme nous l’avons vu en haut, elles 
jouissent des propriétés similaires, de sorte que la puissance r dans la définission (1.54) 
indique une régularité dans la direction de x3 correspondante à celle de H'"(—b, a). 

De manière analogue nous définissons l’espace de Hilbert H°”"(Q*) pour les fonc- 


tions définies sur Q* = R2x]0, a[ par la norme 
CO . 1/2 
Ill sr ço+) = | > l (+12) "a+ +R "IRC, dé) (1.55) 
k=1 
ici A(Ë, k) est la fonction de ë € R? et de k € N\{0} définie par 


_ 2, T 
ü(£,é, k) . Faso (sin Phi U) L2(0,a 


F (0,2) (0) étant la transformée de Fourier par rapport à (x, X2). 
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Système d'équations paraboliques 
linéaires du type : température et densité 


de vapeur avec l'effet de l’évaporation 


2.1 Introduction 


4] Mans le présent chapitre nous allons étudier un système d'équations parabo- 


Su : TRE : PE : D »2 ss 1 
2)  liques linéaires couplées par une condition particulière. L'étude est motivée 


par le phénomène d’évaporation de l’eau de la surface de l’eau liquide, qui intervenant 


dans un grand nombre de phénomènes physiques, joue souvent un rôle important à 
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la grande échelle comme l’évaporation de l’eau de la surface de l’océan (voir par exemple 
[8], [21]) ainsi qu’à l'échelle plus petite ; dans des cas particuliers, des physiciens et des 
chimistes ont proposé des modèles mathématiques bien articulés des comportements 
locaux de l’évaporation de l’eau (voir par exemple [10], [33], [35], [43] ; sur la physique 
de l’évaporation en général, voir aussi [20]). Mais il nous semble que l'étude mathé- 
matiques des équations qui décrivent l’évaporation de l’eau et ses effets n’est pas suf- 
fisamment développée. 

Rappelons les caractéristiques particulières de l’évaporation. Si S est l’intersurface 
entre l’eau liquide et l'air, l’évaporation de l’eau qui se produit sur la surface S est ac- 
compagnée par l'absorption de la chaleur, dite chaleur latente de la transition de phase 
de H20, de sorte qu’il y a une source (négative) de la chaleur concentrée sur la surface 
S. D'autre part, si on admet la présence de la diffusion de la vapeur d’eau dans l'air, la 
quantité de l’évaporation devra être proportionnelle à la composante normale à l’in- 
tersurface du gradient de la densité de vapeur, tandis que sur S la densité de vapeur 
doit être celle de la vapeur saturée. Quant à la diffusion de la chaleur, son coefficient 
ainsi que la chaleur spécifique sont différents dans l’eau liquide et dans l’air. Ces re- 
lations physiques, traduites dans les conditions sur les équations de diffusion pour la 
chaleur et la densité de la vapeur d’eau, constituent les conditions particulières qui ne 
facilitent pas la résolution des équations. 

Dans ce présent travail on considère une équation parabolique linéaire pour la 
fonction inconnue T (température) dans le domaine Q={xe R°|-b<x;<a}etune 
équation parabolique linéaire pour la fonction inconnue IT (densité de vapeur) dans le 
domaine Q* ={xe R°|0 < x3 < a}. On suppose que dans l'équation pour T la “source 


de la chaleur” est concentrée sur {x3 = 0} et proportionnelle à = [r1| II| ss) avec 


X3=E] E 
un €] > 0, ce qui est analogue à l’approximation souvent utilisée dans les travaux pra- 
tiques (voir par exemple [17]) de la loi de Dalton (voir par exemple [1], [37]). On sup- 
pose aussi que la fonction II doit satisfaire à la condition 

Il| = To +] T| 


X3=0 X3=0 
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avec deux constantes %o et &1 ; cette condition estunelinéarisation de la relation II | de — 
T,$(T) avec la densité de la vapeur saturée x,,(T) à la température T. 

Les équations que nous allons considérer sont linéaires et les méthodes de la ré- 
solution des équations de ce type sont bien établies (voir par exemple [13], [23], [2], 
[30]). Mais notre problème exige une élaboration non indifférente à cause du couplage 
particulier mentionné ci-dessus et aussi à cause de la discontinuité du coefficient pour 
l'équation de la chaleur. Les aspects généraux de cette dernière problématique ont été 
investigués par plusieurs auteurs (voir par exemple [22], [4] [5]), mais nous allons uti- 


liser une autre méthode directement utilisable pour la résolution de notre problème. 


2.2 Système d'équations à proposer 


Nous précisons d’abord le domaine dans lequel nous considérons notre problème. 
On pose 


Q = {(x1,22, x3) € R°| —- b< x3 < a}, (2.1) 


où a et b sont deux nombres strictement positifs. Nous définissons aussi la partie infé- 


rieure Q= et la partie supérieure Q* de Q, 
Q = {(x1,%2,x3) € Q[x3 <0},  QŸ={(x,%2,%3) € Q]|x3 > 0}. (2.2) 
Dans le domaine Q nous considérons l'équation 
CopOrT = V-(KVT) + wô(x3), (2.3) 
tandis que dans le domaine Q* on considère l'équation 


Ô,11= yoAIl. (2.4) 
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Dans les équations (2.3)-(2.4), yo est une constante strictement positive, tandis que 


Cye = eo KEK dsO0, Ge ce x=x® dans Q”, (2.5) 


ch, «®, cl, 2) étant des constantes strictement positives; le symbole 6(x3) dans 


(2.3) désigne la delta de Dirac par rapport à x3. Pour y nous supposons qu’elle est une 
fonction définie sur R? et doit satisfaire à l'équation 
Nice, = |, -0 


W=Yi ; (2.6) 


€1 


où y. et €, sont des constantes strictement positives satisfaisant à la condition 0 < €, < 
a. 


Pour les fonctions inconnues T et II nous posons les conditions aux limites 


ThsieTs The Ts (2.7) 
IT |x3=0= To + X1T |x3=0 (2.8) 
D lxs=a= [a (2.9) 
et les conditions initiales 
T |:=0= To(x) x € O, (2.10) 
I |:=0= Ho(x) xeQ*. (2.11) 


Dans (2.7)-(2.9), T_p, Ta La %o0 et di sont des constantes. Pour Ty(x) et Io(x) nous 


supposons les conditions de compatibilité 


Tot-b)=T_p  To(a)=Ta 


Io(0) = %o + &1 To(0),  Ilo(æ)= Il, 
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2.3 Une variante de série de Fourier 


Pour résoudre l'équation (2.3) nous aurons besoin d’une version particulière de 
la série de Fourier, que nous présentons ici dans une forme qui ne fait pas référence 
spécifique à l’équation (2.3). 

Soient L, et L: deux nombres réels strictement positifs (dans l’application à notre 


problème on aura L; = bet L, = a). On considère le domaine en une dimension 


I={xEeR|-L,<x< Lo}. 


On considère deux fonctions strictement positives v(x) et x(x) constantes dans chacun 


des sous-intervelles 7, =] -— Li,0[ et 2 =]0, L[ 
v(x)=v1 sixel, V(x)=v sixeb, (2.12) 
K(X)=Kk1 Si XEl, Kk(xX)=k2 sixeb (2.13) 


avec des constantes strictement positives V1, V2, K1, K2 (on ne suppose pas de condition 
NS 2 2 V] V2 . _ # = à 
analogue à (3.1), donc en général + # +). Au point x = 0, on peut choisir quelconques 


valeurs de v(0) et k(0), donc on peut prendre par exemple 


Vi+V2 


vo= EE, DE 


2 » 


ce qui n'influence pas le résultat. 
Nous allons d’abord construire une série de Fourier relative aux deux espaces de 
Hilbert 


L£(D ={u:1—R, mesurables | [rcoiucoiax < ©} (2.14) 
I 


et 


H\(D={ueC(TR) | frcotu Got ax < oo, u(—Li) = u(Li) =0} (2.15) 
I 
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munis respectivement du produit scalaire 


(u,V)12(D = frcoutovcoax (uv) np = fxcou'to v'(x)dx, 


et de la norme 


lullz2cn = 1/U, U)r207 ul fin = 1/4, U) fn: 
40) #0) k(D x () 


Définissons pour À > 0 les fonctions BQ), YA), R(A) comme suit : 


_ K2V2 Av: . K] T° 
= — L —, 2.16 
BA) arctg| er | / ra 1)) si0 < 1 < RE 7 (2.16) 
S 1 l 
fD=UM-S)T sil LE h=t2e, 


BU) = (n= >yn + arcte| _ 21g| 5) 


1VI K] 
1 
Ur I <A < + Dr n=L2, 
Avi 
1 cos|1/—L; 1 
FU) = cs(pu) 4m, n=l2…, (2.17 


= Si 
V AK2Vo COS BQ) vil 


7 = 


l 
/AK2V2 ViL 


RQ = BU) + 4) 22. (2.18) 
K2 


La fonction A(À) a la propriété suivante : 
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Lemme 2.3.1. La fonction RQ) définie dans (2.18) est continue et strictement crois- 


sante et vérifie les relations 


lim A(À) =0, lim A(À) = oo. (2.19) 
ÀA—0* À—00 


Démonstration. Il résulte immédiatement de la définition (2.18) de Ah) et de la défi- 


nition (2.16) de B(À). 


Le lemme (2.5.1) nous permet de définir une suite {AK}, par 


Ae=h (kr),  k=1,2,-.….. (2.20) 


En utilisant la suite {Ar} on va définir les fonctions e£(x), k = 1,2,---, qui vont for- 


CO 
k=1 
mer une base orthonormale de l’espace de Hilbert L£(1) ayant des propriétés utiles 


pour résoudre notre problème. 


Proposition 2.3.1. On pose 


ex (x) = AS (2.21) 
où 
__. = sin | AL + vx) NP ll _. 
Fupsin (BA + 2x) S0=r< 1. 
Alors 


(A) {ex} est une base orthonormale de Le (1 ; 
(B) {ek}®, est un système orthogonal complet de HA (I); 
(EC) on a 


D 
lex = Â& 
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(D) 


d d 
rire) lire; Vrel=tDr0NnI0 ET 
dx dx 


etK(X) Ler(x) est continue Sur [—-L:,Lo] ; 


(E) on a 
sup lex(x)l < K, 
—Li<x<Li 
l l V/K1 VX 
Ki = 2max| » | ESS) 
VVIiK1 VV2K2 VL; V Lo 
(F) on a 


Motk- 1) < A4 <Motk+1)}},  Mo= | 


TT | 
ENERENES 


2.4 Démonstration de la proposition 2.3.1 


La démonstration de la proposition 2.3.1 s'articule en plusieurs étapes. Nous pro- 


cédons donc en démontrant d’abord des lemmes et ensuite la proposition. 


Lemme 2.4.1. Il existe un opérateur linéaire A de L?(1) sur H}(1) qui à chaque u € 


LE (1) associe Au € H 1 (D) vérifiant la relation 
V K 


L'opérateur À, considéré comme opérateur de L£(1) dans lui-même, est un opérateur 
linéaire auto-adjoint et compact et ses vecteurs propres êx correspondants aux valeurs 


propres Uk, 


AËk = Lkêk, (2.24) 


forment une base orthonormale {4}? , de L£(T). Les valeurs propres x vérifient 
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la relation 


DU M POUR A2, urx—0 pour k—c. (225) 


Démonstration. Puisqu'on a 
A 
plz <Clepläon VpEe H,(D 


avec une constante C et que donc le produit scalaire (u,@)2(N peut être considéré 
comme fonctionnelle linéaire sur @ € H}(1), d’après le théorème de représentation de 
Riesz, pour chaque u € L2(I) il existe un élément U € H1(I) et un seul qui vérifie la 
relation 


(0,9) fan = Ph VpE HO. (2.26) 


En posant Au = U, on définit l'opérateur À de L2(1) dans H}(1). 
Cela étant, le lemme résulte de la théorie bien connue sur l'opérateur linéaire auto- 


adjoint et compact et de la définition des espaces L2(1) et H}(1) (voir par exemple le 


lemme 1 du $ 1 du chap. IV de [28}). 


Lemme 2.4.2. Soit {&k}®, la base orthonormale de L?(1) définie dans le lemme 2.4.1. 


Alors pour tout k on a 


dd d 
pe(re à) = —Vêk dans L°(1). (2.27) 


Démonstration. Comme A+ = uye4, en posant U = Ay, d’après (2.26) on a 


d mx (Lao) Lptodx = Î va ()p (dx. 
I dx dx I 
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Compte tenu que H}(1) est dense dans L?(1), cette égalité définit la dérivée généralisée 


æ (K(x) LC (x)), qui vérifie (2.27). 


Lemme 2.4.3. Soit {ex} ® , la base orthonormale de L?(1) définie dans le lemme 2.4.1. 
Alors pour tout k la fonction x (x) +-8K(x) est continue et la dérivée à gauche et la déri- 


vée à droite de ë(x) au point x = 0 ont le même signe. 


Démonstration. Il résulte immédiatement du lemme 2.4.2. 


Lemme 2.4.4. Soit {&k}®, la base orthonormale de L?(1) définie dans le lemme 2.4.1. 


Alors pour tout k, e,(x) a la forme 


00 = sa 


= ————, (2.28) 
lykllzzcn 


où yK(x) est la fonction donnée dans (2.22). 


Démonstration. On considère la famille de système d'équations linéaires avec un pa- 


ramètre À > 0 


d 1 
= J(X) = = z(x) 
dx K(X) (2.29) 
Lz(0 = -AV(x) 7H 
avec les conditions initiales 
y(-Li) =0, z(-L;)=l. (2.30) 


Par des calculs élémentaires (même si un peu longs) on constate que, si (À) et ÿ() 


sont les fonctions définies dans (2.16) et dans (2.17), les fonctions y(x) = y(A; x), z(x) = 


z(; x), 
1 sin LUS PRES AE si —-L<xXx<0 
yAx=d Van M Hs) (2.31) 
FU) sin (B() +1/%23) si0<x<L 


Khadidja HALLACI 30 Université 8 Mai 1945-Guelma 


Chapitre 2. Système d'équations paraboliques linéaires du type : 


cos| Pr + 4/2 x) si —-L<x<0 

zx) = : (2.32) 
FA) VAk2V2 cos (Ba) + 2x) si 0<Xx< Lo 

constituent la solution unique du problème de Cauchy (2.29)-(2.30) dans la classe 


C([-L1, L])). En outre, comme on le voit immédiatement, on a 


d d 
5 (KO yG) = Av (2 y (2.33) 


presque partout dans ] — Li, L,[ pour tout À > 0. 


D’après la définition (2.16) de B(A),on a 


K] T° 
Ba, HS al 
PAS NL TE ls 2 
AE 2) JF; ri= 
=] me RE Dr, A+ Sr? [, (2.34) 
vil vil 


ce qui implique que 


Av: 
cos A1) 


_ >0 
cos B(À) 


n=1,2;:::; 


donc, d’après la définition (2.17) de 7(À), on a 7(À) > 0 pour tout À > 0. 


Or, d’après la défintion des fonctions A(À) et y(; x) (voir (2.18), (2.31)ona 


y L) = FO) sin(A (À). 


Donc 


y(O; Li) = 0, siet seulement si À = h (kr), keN\{0}, (2.35) 


c'est-à-dire, si et seulement si À = À4 avec un des nombres 1} définis dans (2.20). 
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En retournant à e4 (x), de (2.27) on déduit que, si on pose 


(k) = ee Y(x) = - ëK(x), Z(x) = - ee 
Cz = Ki 00). , PJ(X)= ——28;(x), Z(x) = KO) 780), 


—Li Cz(K) Cz(K) 


alors les fonctions ÿ(x) et Z(x) satisfont au système d'équations (2.29) et aux conditions 
initiales (2.30). Donc en vertu de l’unicité de la solution du problème de Cauchy (2.29)- 
(2.30), on a 


YQ) = y, x), 2{x) = 2, x), 


d’où, en posant y£(x) = y(Ag, x), on a 
er (x) = cz(K) yr(x). 


Comme lekllz2cn =],ona 
l 


Ca) = ——. 
lYkllz2çcn 


Ainsi on a établi (2.28). 


Lemme 2.4.5. Soit y(A; x) la fonction définie dans (2.31). Soient 4 les nombres défi- 
nis dans (2.20). Alors pour tout k, la fonction y(Ax; x) = y£(x) admet k — 1 zéros dans 


l'intervalle ] — Li, Lf. 


Démonstration. De la forme de la fonction y(; x) donnée dans (2.31) et de la propriété 
de la fonction B(À) (en particulier (2.34)) on déduit que le nombre de zéros dans l’inter- 
valle ] — Li, L[ est une fonction croissante et qu’à chaque point À = 1+ il croît d'unité. 


D'autre part, on a 
ou YG;x)=x+1l pour xe[-Li, Lo]. 
—0+ 


On en déduit l'affirmation du lemme. 
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Maintenant nous démontrons la proposition 2.3.1. 


Démonstration. Du lemme 2.4.4 résulte 


êk = CH 
et de (2.33) et de (2.27) on déduit 
l 
Ag = —. 
Uk 


Cela étant du lemme 2.4.1 découlent les propriétés (A), (B), (OC) et (D). 


yk( 
lykll LD 


Nous rappelons que e;(x) est de la forme eg(x) = avec 


AkVi 
K] 


yk(X) = sin | (Li + ») pour —-L<x<0, 


1 
V AgViki 


= + =: Î Ave L O0<x<L 
Pa Ein | 6 Ua 2) pour U=£ X= Lo. 


: d ge d : d sd A 
Comme Jim ki 75 €x(0) = Jin K2 3 €x(X), Jim ax ex (x) et pus 4Lex(2 ont le même 


signe, ce qui implique l'existence d’un x; € 2,01 tel que [sin(\/ An (Li + x))] = 


Agv2 
K2 


ou d’un x, € [0, L2 | tel que | sin | (X1 + L2))| = 1. Dans le premier cas on a 


Lo X1 
role = vomcoaxz [ nincoitax= 
LI AL] 


1 fa Ja 
: Î sin? | PU Gp à x) dx L 
AK Li K] 


1 É sf Er +2) : 


_ AK] 2 2 2 AgVi —Li à 
K1 


l X1 Li Li 
= —(—+—)> à 
AxKi 2 2 AK] 


De manière analogue dans le deuxième cas on a 


Li 


F) 
IykCO 207 = 11. 
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Comme 


1 
sup [eg(x)|=—max( sup |yk(Dl, sup |yk (dl), 
—L1<x<Lo lYkllz2cn —L1<x<0 O<x< Lo 


on en déduit la propriété (E). 
On rappelle que, d’après le lemme 2.4.5, eg (x) admet k — 1 zéros dans ] — Li, Lo. 
Désignons par gq le nombre de zéros dans ] — L1,0] et par r celui dans ]0, L2[ (de sorte 


que qg+r=k-1).Alorsona 


V 

CES M, <(g+Dx, 
K] 

TT<; ARE <(r+1)7, 
K2 


d’où 
V1 V2 
n&-1)=r(q+r)< Va, + 212) <n(q+r+2)=7(k+1). 
l 2 


On a donc 


É 


Motk— 1)? < A < Motk+1},  Mo= | 


TT 
[v [v 
L; 


La propriété (F) est démontrée. 


2.5 Solution stationnaire 


En revenant au système d'équations (2.3)-(2.4), nous considérons avant tout la so- 
lution stationnaire du système d'équations (2.3)-(2.4) avec les condtions aux limites 


(2.7)-(2.9) et avec y définie dans (2.6). Plus précisément, en écrivant z au lieu de x3, 
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il s’agit de trouver les fonctions T(z) et II(z) qui satisfont aux équations 


& 5 4 

- KT) =Wô(z) pour -b<z<a, (2.36) 

d? 
100 pour 0<z<a (2.37) 

avec y définie par 
I). —-1) 
y = y lee = Mle-o (2.38) 
€1 


et aux conditions aux limites 


Thiers Tete (2.39) 
I {2=0= To + 1T |2=0, (2.40) 
I |220= [y (2.41) 


Les deux équations (2.36) et (2.37) sont couplées par la définition (2.38) de y et la 
condition aux limites (2.40). 
Pour déterminer la solution du problème (2.36)-(2.41), on introduit la fonction 


T,(2) continue définie par 


_ GrT _T 
Top + TETE (2 + b) pour -b<z<0 . 
= be . 
‘ Teri. 0<z< 
a me (4-2)  pour0<z<a 


pour laquelle K AT (2) a la même valeur dans ] — b,0[ et dans ]0, al. 


On définit en outre 6, par l'égalité 


_ A _ 
Oo = LT, - (Go + CT e(0) + Oo)))}, (2.43) 
c'est-à-dire 
— 7140 FH = =T 
00 = ——{(lls- 70 - &1Tc(0)), 
a+Y1A0o@1 
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où 


_ CO 0 2 
Us 
k=1 k 


On a la proposition suivante. 


Proposition 2.5.1. La solution du problème (2.36)-(2.41) est donnée par 


Ta — (To + &1(T (0) + O0)) ex(0) 


= ne 8 = k 2.44 
T(2)= Ti(2)+ 2 Cher(z), =" . À (2.44) 
M ne oo ne REP (2.45) 


a 


Nous précisons que la proposition 2.5.1 ne suppose pas la condition (3.1). 


Démonstration. Il est clair que la solution (formelle) du problème (2.37), (2.40), (2.41) 


est 
Ty — (To + 1 T (0 
Hé Sn a top ee More) (2.46) 
a 
Donc, en substituant cette relation dans (2.38), on obtient 
U|,, lo  I,-(Go+41T(0) 
= PS 2 y TE y (T (0). (2.47) 


€] a 


On remarque que K(2) LT (2) a une valeur identique dans ] — b,0[ et dans ]0, al, 


ce qui nous permet d'écrire 


(kw PT) =0 dans ]-b,al. (2.48) 


Nous posons également 


@(z) = T(2)-T,(2). (2.49) 


Alors le problème (2.36), (2.39) se réduit à 


= )) = Wo(T(0))ô(z) (2.50) 
A de — #0 Fe ‘ 
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O(-b) = O(a) = 0. (2,51) 


En posant 


@(z) = Ÿ_ cPer(z) 
k=1 


et en intégrant les deux membres de (2.50) multipliés par e-(z2) sur [-b, a], en vertu de 


(B), (C), (D) de la proposition 2.3.1, on a 


æ d d 
AkC® = [ac(- (KO) Jdz = Wo(T(0))ez (0). 


C'est-à-dire, la solution du problème (2.50)-(2.51) est donnée par 


00 0 
e= ete,  =potTo) O. (2.52) 
k=1 À 
Donc la valeur de @(0) est 
. . ce 0 2 
EU =ptTO)20, oz > EU, (2.53) 
k=1 A 
De (2.47) et de (2.53) on obtient 
A — e 
@(0) = T (Ts — Go + &1(T(0) + O(0)))). (2.54) 


Donc, d’après la définition (2.43) de @o, on a 


@(0) = @o 


et donc 


T(0) = T,(0) + Oo. 


En substituant cette expression dans (2.47), (2.52), (2.46), on obtient (2.44) et (2.45). 
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système d'équations liniéaire 


) À _VANT d'envisager le problème complet (2.3)-(2.11) dans Q € R, nous allons dé- 


montrer, dans la section 3.1, l'existence et l’unicité de la solution dans le cas où 
le domaine se réduit à l'intervalle ] - b,a[€R, c’est-à-dire, dans le domaine d’une di- 
mension spatiale. Pour le problème complet dans Q € R°, nous démontrons le résultat 


d'existence et unicité de la solution sous l'hypothèse 


D 

ve __ Cve 
= —— (3.1) 
KO 2 


38 
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Il est clair que la condition (3.1) est assez restrictive. Mais, comme dans le cas du do- 
maine d’une dimension spatiale le résultat s'obtient sans cette restriction, nous espé- 
rons obtenir, dans un futur prochain, un résultat analogue sans supposer la condition 


(SL), 


3.1 Cas du domaine d’une dimension spatiale 


Considérons le problème réduit à une dimension spatiale. Il s’agit donc du système 
d'équations 
CopÔtT = 0,(K0,T) + WÜ(z) dans R;x]—b,al, (3.2) 
011=Yo0I1 dans R;x]0,al, (3.3) 
avec 


IT(E1) — IT(0) 
— 


: (3.4) 


les conditions aux limites et initiales sont (2.7)-(2.11), dans lesquelles les fonctions ne 
dépendent pas de (x1,x2). Dans ce problème on ne suppose pas la condition (3.1). 
Dans la section précédente nous avons déjà montré l'existence et l’unicité de la 


solution stationnaire de ce problème, que nous notons (T,4(2),I1,,(2)). Posons 
O= 0,2) = T{t,2)-Tse(z), (3.5) 


n = (6,2) =I(t, 2) —T,4(2) — 19 (6, 2). (3.6) 


Comme 62II,, = 0 (voir (2.37)), de (3.3) découle 
On —Yo05n = -@&1(0:9-7Y0059) pour 0<z< a. 


On remarque que, en vertu de (3.2), (2.5) et (3.5),on a 


_ | +0 ; 
—&1(0;9 — 700,0) = (=; -70)05 pour 0<Zz< 4. 
Coop 
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Ces relations nous permettent de réduire le système d'équations (3.2)-(3.4) à 


Cvodr0 — 0,(K0,0) = qÔ(2) dans ]- b, al, (3.7) 
0nm-7Yod:n=f dans ]0,al, (3.8) 
q = q{) ne, (3.9) 
f= a( — 0029. (3.10) 
Cig 


Les conditions aux limites sont 


O(E, —b) = O(F, a) =0, n(t,0) = n(f, a) =0 Vt2>0, (3.11) 


tandis que les conditions initiales sont 


9(0, 2) = To(z) — T'sr(2) = Po(z) pour -b<z< a, (3.12) 


n(0, z) = [o(z) =I1.,;(2) — &100(2) = 7o(2) pour0<z< a. (3.13) 


Pour les fonctions u définies sur [-—b, a], en utilisant les fonctions e- définies dans 
(2.21)-(2.22) avec L; = b, L, = a, on définit les coefficients de Fourier &(k) et la norme 


par 


ul. 


NL 


A = (uen) ce Mr, =(D AE), 8.14) 
k=1 


où Le (—b, a) est défini comme dans (2.14). D'autre part, pour les fonctions u définies 


sur [0, a], en utilisant les fonctions gz(2) = VE sin(?K2) on définit 


CO 1 
A0 = (geo lle, = (2 ERP). (G.15) 
k=1 
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Proposition 3.1.1. Soit 1 < r < Ÿ. Supposons que 99 € ie etno € H%.,,- Alors, quel 


que soit t > O, il existe une solution (8,1) et une seule du problème (3.7)-(3.13) dans la 
classe 


de L°(0, 640) NL 0, 6H po) 


He O 6 nd) (3.16) 


Démonstration. En utilisantles coefficients de Fourier définis dans (3.14)-(3.15) (même 
si on utilise le même symbole #(k) pour les deux sens différents, il n’est pas difficile de 


les distinguer du contexte), on peut transformer les équations (3.7)-(3.8) en 


TB, R) + A40(4, k) = q(ner(0), (3.17 


d . RE _ 
a 0 +705 k ñ(t, k) = f(E, k). (3.18) 


En multipliant les deux membres de (3.17) par AE, k),ona 
1 Eee 
OA PE, k)) + 4792 (1,0) = A TA (t, D q(P er (0). 


Comme 


leg(0)| < Ki VkEe N\{0} 


(voir (E) de la proposition 2.3.1), on a 
r-14 k 1 ra k 2 K} r—2 2 
A4 d(t, )g(er(0) < ASIE, )| Fo A lg). 
Or, comme r-2<-— D en vertu de (F) de la proposition 2.3.1 on a 


(ee) 
a À = Ka < oo. 
k=1 
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On en déduit que 


SUBI +8,  <KKUqO. 


He 


Quant à |g(#l*,ona 


2 
qui = Éntent = É|S Sn by/Esnkenf < 
ET 


1 


<Ky Le AG, D = Kln(, JIÉ 


(0, 4 


RSA 


Pr 


47,2 
GET k=1 


En adjoignant les deux inégalités, on obtient 


TBE + ON < KÉKAKGIMA (8.19) 


(0,a) 


HE 


D'autre part, si on multiplie les deux membres de (3.18) par k*"-?(#, k) et on fait la 


somme par rapport à k, on obtient 


2 2 
UE 1 yo Æ rl . < KIf(, yes (3.20) 
où 
1 a 
Ko = —— 
Yo 
Si on définit 
D=e "9, G=e “in q=e "ta, f=e"!f (3.21) 
avec un & > 0 (à choisir dans la suite), de (3.7)-(3.8) on obtient 
Crp019 — 0,(K00) + WCypŸ = G6(z) dans ]-b,al, (3.22) 
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Ô-Y002+wÿ=f dans ]O0,al. (3.23) 


En procédant d’une manière tout analogue à l’obtention de (3.19)-(3.20), de (3.22)- 


(3.23) on obtient 


d = 2 = 
1912, +10lr  +2wl012,, <K2KaK IAE , 3.24 
ai) la, I Ur I PE 1 KA Ur (3.24) 

d ns _. Se 

— Fi . +20|ñ1|2,1 <K no. 3.25 
aire 072 Ir, Ag 21 gr (3.25) 

Soit 9 une fonction appartenant à la classe 
Helen niet) (3.26) 


Alors, en substituant ® à la place de ® dans la définition (3.10) de f, on définit f, qui ap- 
partiendra à L(0, f; H . En substituant f à la place de f dans (3.8), on peut résoudre 


l'équation (3.8) (avec f) et on obtient 
neo EH nr (0, 5H) 


(voir (3.20)). Ensuite, en définissant g par (3.9), on résout l'équation (3.7) et on obtient 
0, qui, en vertu de (3.19), apprtient à la classe Y2. C'est-à-dire, on a définit un opérateur 
Gi qui, à 9 € Yo, associe la solution 9 € Y) de l'équation (3.7). 
Si on définit f = e *{ f et si on rappelle la définition de f, on voit aisément qu’il 
existe une constante K telle que 
fl <Krlolr 3.27 
MAS nl [Fa (3.27) 
Considérons maintenant deux fonctions À et % appartenant à la classe Yo. En 


définissant 01, 0», 01, 02, 11, 2 de la manière indiquée ci-dessus, on pose 


Ô=9,-#%h, À =T1 -2, O = 91 — 0. 
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Comme les équations (3.22)-(3.23) sont linéaires, de la même manière que (3.24)- 


(3.25) et (3.27) on obtient 


TIGE +Oly +206, <KIKAKGI AE 1, (3.28) 
(-b,a) (-b,a) Hé a 
T° 
MA 1 +70 Age a + 20 Alryr < KKÇ NO l4y . (3.29) 


Comme 4?! < 64? + 422 et donc | HI? , < 61H12 7 + 5H, r1 Pour tout 
F3 0,4) 


(0,a) 
ô > 0, en multipliant les deux membres de (3.28) par une constante À à a et 


en faisant la somme de ces deux inégalités, on obtient 


d ee 
ne ar (ME s + Are 1) + ANG + 
ia, 
+ ll, +240: _ +201 Hlr 1 
< OAKŸKaKql HG + LREK lÉI2, + KKrlOI2 (3.30) 
_ Ho HGaw 46 1 #4 Fe 7 ‘ 
Si on choisit À =2K2Kf, Ô = RER 0 = ÀK?KaKa, l'inégalité (3.30) se réduit à 
d ” = 
AI® +1 + AIOlÉr  +Yo—1lAlér + 
dr — (AI lys I Vy- 1) I Ur 052 I Ur, 
À = 
+2A0|6l + ol Ale _ < 516 (72 . (3.31) 
L'inégalité (3.31) implique que l’approximation successive On: n =1,2,.., pour 


le problème (3.22)-(3.23) (avec les conditions initiales (3.12)-(3.13) pour tout n avec 


00 € H7° . et no € Hi, a) converge dans l’espace de Banach Ÿ9 muni de la norme 


1/2 
IO1g, = Ce GHT., : a HER 6H 4 a) +20 012 GA, : ” 


Comme les équations sont linéaires, la limite 9 de l’approximation successive d() Sera 


la solution du problème (3.22)-(3.23) (et, par sa construction, donnera aussi 7). 
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En posant 9 = e”! ô, n = e*!, on voit aisément que (5, n) sera la solution du problème 


(3.7)-(3.13). Leur appartenance aux classes (3.16) résulte des estimations établies ci- 


dessus. 


La solution (9, n) du problème (3.7)-(3.13) étant trouvée, par les égalités (3.5)-(3.6) 
on peut construire les fonctions T(f, z) et II(£, z) qui satisfont aux équations (3.2)-(3.4) 


et aux conditions (2.7)-(2.11). 


3.2 Cas du domaine Q CR° 


3.2.1 Transformation des équations 


Dans le cas où le domaine en considération est Q © R* défini dans (2.1), nous allons 
démontrer notre résultat sous hypothèse (3.1). Dans la démonstration, même si nous 
suivons l’idée générale adoptée dans le cas du domaine d’une dimension spatiale, la 
structure des équations aux dérivées partielles dans un domaine de dimension 3 exige 
une élaboration technique non indifférente. 


De manière analogue à (3.5)-(3.6), nous posons 
9 = 0(t,x) = T(L,x) Ts), (3.32) 


n = (6,2) = T6, x) — Ts (x3) — &1O(E, x). (3.33) 


Comme AIl,; = 0 (voir (2.37)), de (2.4) s'obtient 
05n—YoAn =-«1(0;9-70A9) dans Q*. 


Or, de (2.3), de (2.5) et de (3.32) on déduit que 


L. : x) 
—%:1(0,9— 7oAË) = (> — yo) A® dans Q*. 
C6 
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Ainsi on peut réduire le problème (2.3)-(2.9) au système d'équations 


Cvpdr0 — V-KVD = qô(x3) dans Q, (3.34) 


0;n-YoAn=f  dansOQ*, (3.35) 


avec les conditions aux limites 


Ÿ |x3=-Dp= Ÿ lxy=a= 0, I] lx3=0= 1 lx3=a= 0, (3.36) 


et les conditions initiales 


ne) |5=0= To = La = 0 dans Q, 1 |=0= IL — 11: — &100 = 10 dans Q*, (3.37) 


où 
(£, X1, X2, €1) 
q= qu, n,0) = NT 2 do (3.38) 
1 
_ RO 
De a(— —yo]A®. (3.39) 
C 
vo 
On rappelle que, comme nl x=0 = 0 0na 
L Nlxz=e : Nlxz=e —1]|x;=0 
q=Yi1 el = 1 e , 


ce qui fait correspondre (3.38) à (2.6). 

Pour obtenir des estimations utiles à la démonstration de l’existence et de l’uni- 
cité de la solution, nous allons utiliser la transformée de Fourier en (x1,x2) et le dé- 
veloppement en série de Fourier par rapport à x3. En rappelant (3.1), nous réécrivons 


l'équation (3.34) dans la forme 


ee (3.40) 
MTV Ge 6x5 6x | 7 9” 


Khadidja HALLACI 46 Université 8 Mai 1945-Guelma 


Chapitre 3. Existence et unicité de la solution 


où 
: au) : KP) 
7 D 2" 
vo Cp 


Comme en vertu de la proposition 2.3.1 on a 


-[ 3 [x x )) = Ar(er, 9) 
(24 3 Ôxz Ôx3 X3 — Ak ex; Li (ba), 


en faisant le produit scalaire des deux membres de l'équation (3.40) avec e-(:-), de leur 


transformée de Fourier par rapport à (x, 2) € R? on obtient 
0 (4,Ë, k) + (hlél? + Ak)O(r,é, © = Gt, der (0), (3.41) 
où = (1,62), IF = 85 +65 et 
LED = Fox) Gr 92 Car AO F(rux) 


F(0,x) () étant la transformée de Fourier par rapport à (x1, x2). 
D'autre part, en exprimant les deux membres de l'équation (3.35) en la transformée 
de Fourier par rapport à (x1, x2) et en le développement en série de Fourier par rapport 


à x3 sur la base ésinEkene on a 


2 
GE, + ro[le + LE )AtE, 0 = F6. (42 


Ici et dans la suite, pour ne pas alourdir les notations, nous désignons par le même 
symbole ” la transformée de Fourier par rapport à x. et x2 ainsi que le développement 
en série de Fourier par rapport à x3 d’une fonction sur la base {e}®, et sur la base 
(2 sin(Fkx3)}® .. Même si le symbole © désigne différentes transformations, il ne 
sera pas difficile de comprendre du contexte ce qu’il signifie. 


Pour démontrer l'existence et l’unicité de la solution du problème (3.34)-(3.37), 
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nous posons 


d=e "9, G=e "mm ÿG=e"’qg f=e"*!f (3.43) 


avec une constante w > 0, que nous allons choisir dans la suite. Alors on a 


” Ds Das. ds 0e = — 
Copbb— Rovo(s 7 + ou an (ax ©) + 06000 = qô(x3) dans Q, (3.44) 
0ÿ-7yoAñ+wÿ=f  dansO*. (3.45) 


En outre, de manière analogue à (3.41)-(3.42) on a 
d,Ô(,E, 1 + (REP + Ax)O(, 8, + w(L,€, 0 = 4(t, Der), (3.46) 
SZ 2 où 2\= S = 
o,E,0 + ro[iel + LR JF(,E, 0 +oû(,8,0 = (8,0. 47 


3.2.2 Espaces fonctionnels et estimations 


Nous introduisons les espaces fonctionnels que nous allons utiliser pour montrer 
l'existence et l’unicité de la solution de notre problème. 

Pour les fonctions définies sur R? nous utilisons les espaces de Sobolev usuels H°(R?) 
munis de la norme 


1/2 
lullyrque = | [ A+ AG PdE) | (3.48) 


Pour les fonctions définies sur Q* = R?x]0, a[ nous définissons les espaces H°”(Q*) 


par la norme 
S 2\s-r 2, p2ris 20 
ares (À [a+ arr éiac ae) 649 
k=1 


tandis que pour les fonctions définies sur Q = R? x ]—b, a[ nous définissons les espaces 
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H°”(Q) par la norme 
2\S—Fr 2 TA 2 1/2 
ll ro = (E J OH) "(++ a) RG HP dé) , 6.50 


où 14 sont les nombres définis dans (2.20) (avec L; = b et L) = a). 


Cela étant, nous considérons deux nombres réels r et s satisfaisant aux conditions 
3 
1<r< 7’ S>T. (3.51) 


On va utiliser les espaces H°-2(R2), H°-l"-1(Q), etc.., qui sont définis selon (3.48)- 
(3.50) (avec la substitution de 5 — 3 au lieu du nombre générique s, etc...). La motiva- 
tion de l'introduction des espaces H°'(Q*) et H°”(Q) avec s > r est que pour l’inter- 
prétation physique il est important d’avoir la bornitude de la restriction à {x3 = 0} de la 
température T, qui sera obtenue en supposant que 3 > S. 

Dans la suite nous allons établir des estimations de la solution 9 de l'équation (3.34) 
(ou (3.40)) et de la solution 7 de l’équation (3.35), estimations en fonction de g et de f 
considérées à titre provisoire comme données. Puis, on établira aussi des estimations 
de gq et de f définies par (3.38) et (3.39). Pour leur démonstration, nous allons utiliser 
l’idée du théorème de trace illustrée dans [26]. Toutefois, comme nous utilisons des 
espaces de Sobolev non usuels, nous démontrons les inégalités d’une manière directe 


sans recourir aux résultats connus des théorèmes de trace. 


Lemme 3.2.1. Soit Ÿ une fonction satisfaisant à l'équation (3.40) et à la première 
condition de (3.36). Soient à et q les fonctions définies dans (3.43). On suppose que 
q € L2(0, t;: H°-2(R?)). Alors on a 
ETTE alle à cal. , (3.52) 
Fret) (Q) He 


ne (R2) 


TBE. 1,r— 1(Q) “F cle +20| 817. dr 1(Q) = < Gale. Re) (3.53) 
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où 
À e K2 22-T 
— nie ñ Ci = 2 (1 . 
1+4 Ci (3-2r)V Mo 


À (3.54) 
Démonstration. On multiplie les deux membres de (3.41) par (1 + é|?)°"(1 + hjël? + 
+Ap)" LA (E, 6, k), de sorte que l’on obtient 

l : ne 

CO D CO TE 7° CT 0) 


++) TG + RE + 0 HE + A2) O7 (6, €, k) = 
= (+) "A+ hé +19" (4, E, Dr, Der(0). (3.55) 


Comme À; > 0 (voir (2.20)), on a 
c1(1+ hlël? +) < hlél? + Ag. (3.56) 
D'autre part, d’après (E) de la proposition 2.3.1 on a 
ler(0)| < K VkeN\{0}. 


On a donc 


A+)" + hé +10) TO (L EE, RG, Der (0) < 
< Su +)" + Réf + A) 10 E, R)I2+ 
K2 
Fe HET CLEA + M IG OÉ. 


Commer-2<0,ona 


A+hleP+2p)" 2 <227 0/14 hé + VAT. 
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Donc, en vertu de (F) de la proposition 2.3.1 on a 


œ ; , 2 2-7 © Vl+hléZ 4 
HU A 7 UNS de US CE Te 
L ‘ Ge) L V Mo 
on aenoutre 
CO 2 
et +774 < 
&1 VMo 
Le] 2 2,23 
<[* 1+hl cn 1 = jo 
0 VA 3-2r\ M 
et 
A+hé2 +) 2 < (+ hé". 
On a donc 
CO 2—r L 
LL RER TZ (3.57) 
_ D (+7) 


En intégrant par rapport à 6 = (61,62) les deux membres de (3.55) et en faisant leur 
somme par rapport à k (voir aussi les définitions (3.48), (3.50)), compte tenu des rela- 


tions (3.56)-(3.57), on obtient (3.52). 


L'inégalité (3.53) s'obtient d’une manière tout analogue (voir (3.46)). 


Lemme 3.2.2. Soit n une fonction satisfaisant à l'équation (3.35) et à la seconde 
condition de (3.36). Soient ñ et f les fonctions définies dans (3.13). On suppose que 


fe L2(0,r; H$-2"-2(Q*)). Alors il existe une constante C3 > 0 telle que 
SAT - < GI flÉ 3.58 
7 sr (04) + VoC21nl rs ça+) = 31 ll7rs-2r-2ç0+r (3.58) 


d ZI 2 ZIIZ ZIIZ 7 1 FIL 
dt (M 3s-17-1 ç0+) + VoC21 Frs a+) 76 201 s-1r-110+) < Call ls-2r-2ç0+) (3.59) 


(une expression explicite de C; est donnée à la fin de la démonstration). 
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Démonstration. On multiplie les deux membres de (3.42) par (1 + [é[?)°"(1 + |ë? + 


+k2) (1, ë, k), de sorte que l’on obtient 


l 
3 0(Q + TA + lélÉ + RE (6, , ©))+ 


2 
+yo(L + El) Q + EE + KDE + RU, ë,k) = 
= (1+1é 2) a+ + 2716, €, D f(L,E, 0. (3.60) 


On rappelle les inégalités 
2 


2 

TT 
Co(1+lËél+Kk2) <lé+—Kk2, c=min(1, —), 
2(1 +6] )<{] 2 ( 7) 


a+ +2) 19, ER ft, E,k) < 
0 4 1e + GE 2 + + 12 + LIFE D. 
2 2C2Y0 


Si on substitue ces inégalités dans l'égalité (3.60) et en fait l'intégrale par rapport à 
& = (61,62) et la somme par rapport à k, en raisonnant de manière analogue à la dé- 


_ l 
monstration du lemme 3.2.1, on obtient (3.58) avec C3 = —. 
C2Yo 


L'inégalité (3.59) s'obtient d’une manière tout analogue (voir (3.47)). 


Lemme 3.2.3. Soient q, q et ñ les fonctions définies dans (3.38) et (3.43). Alors, quel 


que soit Ô > 0, il existe une constante C5 telle qu'on ait 


1141 Ë7s-1 qe) = nee + C7) (3.61) 


= ZIlZ 2 
IG 5-1 2 = OT As ço+) + Cols 0) (3.62) 


Démonstration. Comme 


2 2 k 
na en = D GER Esin( Le),  E= (né), 
k=1 4 1 
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on a 


2 sil 2 2 [ Al - 
for msenrars [a+ (pme ei) ar 


On remarque que 


OO OO 1 (ee) 1 
AC, &, k)| A++ KT 219 (t, ë, D? 
(2 Ë <(Z de 2 


et que l’on a 
CO ce) r-1 
D À a 
fi 0 1 HI + KT 
00 7-2 1-2 
0 (V1+lé/2+x)271 (2r—-2)(1+|é/2)71 
On a donc 


[ A+ AU, Een dE < 


2-21? 2\s—r 2,r- 
= A2 a LL DE Da+kr + KT 21 (5, , RP dé < 


2 s—r 2 
sf, 0 + D a+kr + VIA DPdE) + 


«(far Eu+t +) GE BP4E) 


Par conséquent on a 


2 


2 Yi > 
1915-21 2 L _ InC, 5° €1) Îgs-1 pe) = 
l 


= sf,a+ BTE A++ KO IAE, HF dé+ 
k=1 


Ge [ He À + lé + k27-190,6, 0 dE, 
k=1 


où 
2r+1,,4 
2 vs 


Co = — 
À 46 a?(2r —-2)?€1 


On en déduit (3.61). 
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L'inégalité (3.62) s'obtient de la même manière. 


L'estimation de f définie par (3.39) que nous allons utiliser résultera immédiate- 


ment de (3.39) et sera mentionnée directement dans la section suivante. 


3.2.3 Existence et unicité de la solution 


Comme nous avons transformé le problème (2.3)-(2.9) pour (T,ID en le problème 


(3.34)-(3.37) pour (9, 7), le problème se réduit à démontrer la proposition suivante. 


Proposition 3.2.1. Supposons que 90 € H°7l""1(Q) et no € H°7l"-1(Q*). Alors, quel 
que soit t > O, il existe une solution (8,1) et une seule du problème (3.34)-{3.37) dans 
la classe 


de L®(0, r;H° b-1(0)) n L(0, 6 H°'(0)), 


ne L®(0, 6; H°7 LT 1(0Ÿ)) n L2(0, 5 H°7(Q*). (3.63) 


Démonstration. Soit 9 une fonction appartenant à la classe 


Y9 = L®°(0, #;: H°-L"-1(Q)) n L?(0, t; H°"(Q)). (3.64) 
Alors, en définissant 
_  _x@ ; 
f= (— = yo)A® (dans Q*), (3.65) 
Cp 


on a 


fe, H°272(Q*), 


ce qui nous permet de trouver la solution 7 de l'équation (3.35) avec f = f, solution 


qui, d’après le lemme 3.2.2, apprtient à la classe 


ne L®(0, 5; H°° LT 1(QŸ)) n L?(0, r; H°7(Q*)). 
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Si on définit q par (3.38), d’après le lemme 3.2.3,on a 
qe L'O,5 H°R?)) € L2(0,1; H°-2(R?), 


ce qui nous permet de trouver la solution 9 de l'équation (3.34), solution qui, d’après 
le lemme 3.2.1, apprtient à la classe Yo définie dans (3.64). 

Ainsi, cette chaîne de résolution d'équations définit un opérateur G qui, à 9 € Yp, 
associe la solution Ÿ € Y9 de l'équation (3.34). 

Pour trouver le point fixe de l'opérateur G, nous passons aux fonctions ë, ñ et q 


définies dans (3.43). Définissons aussi 


Il 
5 
E 
D 


: 


c'est-à-dire, f définie dans (3.43) avec f = f. On remarque que de la définition (3.65) 


de f et de celle de f il résulte immédiatement que 
IF lgs-2r-2çor) < CO Grey (3.66) 
avec une constante Cf. D'autre part, des inégalités (3.53), (3.59), (3.62) on déduit que 


d 9: TP 12 
PAIE ET + Pl Ésr co) + AC OT Frs co + 


SIIZ 112 SIZ 
+C2Y0l 7 çQ+) + 2A& | 0155-15-11) + 2019 17-104) £ 
+ = +. ee 
= AC ÔIT sr a+) + ACC 75-17-21 (04) + Call fllas-2r-2ç0+) (3.67) 


où À est une constante à déterminer (ô et w sont aussi à choisir). Si on choisit À = 


3%, ô = ei w = AC1Cs, de (3.66) et de (3.67) on déduit 


d ET RIZ à112 
dt (AIDE + lésion) + AC; ON Frs co) + 
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C2Ÿ0 


+ Mlsros + 210181 


2 7 SII2 
HS-Lr-1(0) + AC: Co 8-12 ç0+) < 


AC, = 
= lo (3.68) 


Si on considère deux fonctions 9 et 92 appartenant à la classe Ys, alors en définis- 


sant 1, 02, 01, 02, 11, 2 de la même manière et en posant 


O=d-0, H=mi-ÿe  O=01-0, 
les fonctions 6, H et @ vérifient l'inégalité analogue à (3.67), d’où on obtient 


AIË 2 - + Aa IE, - < Sp, (3.69) 
L®(0,t;,HS-Lr-1(Q)) l L2(0,6H7(Q)) — 3 L?2(0,r,H°7(Q))° : 


L'inégalité (3.69) implique que l'opérateur G qui, à 9, associe Ÿ est une contraction 


dans l’espace de Banach Ÿ9 muni de la norme 


_ _ _ 1/2 
_— 2 2 
Ir, = (AE one . AGE ur con) ‘ 


On en déduit l’existence et l’unicité de l'élément à € Ÿs tel que G(9) = à. A partir de 


cet élément Ÿ € Ÿ) on peut construire la solution (9,1) du problème (3.34)-(3.37). La 


proposition est démontrée. 


La proposition 3.2.1 étant démontrée, par les égalités (3.32)-(3.33) on peut définir 


(T, ID, qui sera la solution unique du problème (2.3)-(2.11) sous les conditions 


To—Tsre HT (O), 


Ho Tlse-@ilTo- Tel], € H° 7). 


Il est clair que la régularité de T et II sera caractérisée par 


T-TeL®(0, t; HP 1(0)) n L?(0, t; H°” (Q)), 
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N-T-miT-Ti|,€L°@,5H 7" (Q)nL (0, 5 H°"7(09) 


pour tout f > 0. 
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Équation de la chaleur avec l’effet de 
l'évaporation pour une gouttelette 


sphérique 


4.1 Introduction 


4] ANS les chapitres précédents (voir aussi [18]), nous avons considéré le système 


C2) de deux équations paraboliques linéaires, modélisant le comportement es- 
sentiel de la température et de la densité de vapeur dans l'air et dans l’eau avec l'effet 


de l’évaporation qui se produit sur la surface de l’eau, et nous avons démontré 
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Chapitre 4. Équation de la chaleur avec l'effet de l’évaporation … 


l'existence et l’unicité de la solution de ce système dans un domaine borné délimité 
par deux plan horizontaux. Pour ce faire nous avons utilisé une variante de la série de 
Fourier appropriée à un intervalle ] — b, a[ avec des coefficients différents dans deux 
sous-intervalles. 

Dans le présent chapitre, en utilisant l’idée développée dans les chapitres précé- 
dents, nous allons étudier le problème analogue dans un domaine sphérique. Rap- 
pelons que l’évaporation de la vapeur à partir d’une gouttelette (pas nécéssairement 
de H20) a intéressé en particulier beaucoup d'ingénieurs dans les applications indus- 
trielles. La recherche théorique dans le cadre de la physique mathématique autour de 
ce phénomène a été, elle aussi, développée par plusieurs chercheurs (voir [7], [34], [32], 
[9], [36], …). Dans ces travaux, généralement, ils considèrent la gouttelette comme une 
sphère. 

Suggérés par ces investigation, nous allons considérer une sphère occupée par l’eau 
et une région autour de cette sphère jusqu’à certaine distance occupée par l’air. Dans 
ce domaine, sous l'hypothèse de la symétrie sphérique, nous allons étudier le système 
d'équations paraboliques linéaires représentant le comportement de la température 
et de la densité de vapeur avec l’effet de l’évaporation, en utilisant une variante de la 
série de Fourier relative au domaine sphérique, construite suivant une idée analogue à 
celle de la série de Fourier introduite dans les chapitres précédents. 


Plus précisément, on considère une gouttelette sphérique occupant la région 


Boa={xER|0<r = \/xt+x8 +8 < a) 


et sa région environnante 


Bab=ixeRla<r=:/x?2+x2+x2 <b}, 
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où a et b sont deux constantes telles que 0 < a < b< co. On pose aussi 


Bo,b = BaU Bab = {XERI0<r=4/x?+x2+ x < b}. 


Dans la région B,p occupée par l’air et l’eau, nous allons considérer l'équation para- 
bolique linéaire (4.1) pour la température T, tandis que dans la région B,,, occupé par 
l'air, nous allons considérés l'équation parabolique linéaire (4.2) pour la densité de va- 
peur II. Pour les conditions de couplage nous allons supposer des conditions (4.3) et 
(4.7) analogues aux conditions de couplage représenter dans le chapitre 2. Nous allons 
démontrer l’existence et l’unicité de la solution du système (4.1)-(4.2) avec la condi- 
tion du couplage et d’autres conditions aux limites, que nous allons préciser dans la 


suite. 


4.2 Formulation du problème à symétrie sphérique 


Dans les domaines 


Qo,p ={xER°I0<r=4/x2+x2+3x5 < b} 
Oap={xeRla<r=\/x?+x2+x2 <b}, 


on considère les fonctions T et II, qui représenteraient la température et la densité de 


la vapeur et qui doivent vérifier le système d'équations 


CopOrT = V-(KVT)+wô(r — a) dans Q,,p, (4.1) 
0/11 = yoAI dans Q,,p, (4.2) 
avec 
y = 7, Tate Mie (4.3) 
€] 
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où Yo, 1 et €1 sont des constantes strictement positives satisfaisant à la condition 0 < 
€, <b—a et ô(r’) désigne la delta de Dirac par rapport à r’ (nous y prenons r’ = r — a). 


On suppose que 


(1) 


Gy=6 pourO<r<a& Cy=c) pour a<r<b, (4.4) 
x = x pour 0O<r< a, x = Kx0) pour a<r<b, (4.5) 
ae. KO), ee, k@) étant des constantes strictement positives. Pour les fonctions incon- 


nues T'et II nous posons les conditions aux limites 


T7, (4.6) 
[r=a= To + MT |r=a (4.7) 
W-2= 1; (4.8) 


Rappelons l'expression du laplacien appliqué aux fonctions à symétrie sphérique 


(dans RS). 


Lemme 4.2.1. Si u € C?(Qo,a) ou € C?(Q4,p) et u ne dépend que de r = \/ x? + x? + x£, 


on a 


Sr ÉN 
auG= TE = u(r), (4.9) 


Démonstration. Comme il est bien connu, il résulte immédiatement des calculs expli- 


cites. 


Donc, sous hypothèse que T(f,-) et II(f,-) sont à symétrie sphérique, les équations 


(4.1)-(4.2) se réduisent à 


Cop0:T =6,(r"k(r)6,;T)+yô(r-a) dans R;*]0,b[, (4.10) 
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011=y00-(r 0,11 dans R;*x]a,bl. (4.11) 


La définition (4.3) de y et les conditions aux limites (4.6)-(4.8), ayant l'expression rela- 
tive seulement à r, restent invariantes. 
Cela étant, dans le présent travail nous considérons le système d'équations (4.10)- 


(4.11), (4.3) avec les conditions aux limites (4.6)-(4.8) et les conditions initiales 
T |;:=0= To(r) pour0O<r< b, (4.12) 


IT |;-0= IHo(r) poura<r<b. (4.13) 


On suppose les conditions de compatibilité 


To(b)= Ty, Hot) =70+@1To(a),  Ho(b) =. (4.14) 


4.3 Espaces L° et H!' avec des poids 


Pour étudier le problème (4.10)-(4.11), (4.3), (4.6)-(4.8), (4.12)-(4.13), nous allons 


utiliser les espaces de Hilbert 


b 
L?(0,b) ={p:(0,b) —R, mesurables | | v(r)r{lg(r)dr < oo}, (4.15) 
0 
> _ b d 
H}(0,b) = {pe L?(0,b) If KO par < co, p(b) = 0}, (4.16) 
0 
. b 
L°(a, b) = {@ : (a, b) —R, mesurables il r*|g(r)l?dr <co!}, (4.17) 
1 . 2 ? 2 d 2 _ : 
H'(a,b)={peL"(a,b)|| r brel 10] dr <oo, @(a) = @(b) = 0} (4.18) 


munis respectivement du produit scalaire 


b b 
(UV), = [ v(rn)r£u(r)v(r)dr, ORDENS = k(r)r£u/(r)v'(r)dr, (4.19) 


Khadidja HALLACI 62 Université 8 Mai 1945-Guelma 


Chapitre 4. Équation de la chaleur avec l'effet de l’évaporation 


b b 
= | rfu(r)u(r)dr, & ça = | rfu/(r)v'(r)dr (4.20) 


et munis de la norme correspondante à ces produits scalaires. 


Rappelons d’abord les propriétés élémentaires de ces espaces de Hilbert. 
Lemme 4.3.1. L'injection de H}(0, b) dans L2(0, b) est compacte. 


Démonstration. Si on définit l'application / 
J(p)G9 = px),  xEeQo,p (ou EQ,,p), (4.21) 


alors comme on constate par des calculs, quelle que soit la fonction mesurable  défi- 


nie sur ]0,b[,ona 


b 
[ JM @U(p) dx = 47 À v(nr£lo(r)édr, (4.22) 
0,b 


b d 2 
Î TGV TG dx = 47 | x(Nr2| = p(n)| dr. (4.23) 
Qo,b 0 dr 


Soit M un ensemble borné de H}(0, b). De (4.23), joint à l'inégalité 0 < k1 < J(K)(x) < 


K2 < ©, il résulte que l’ensemble 
JM = {v = J(u)lu e M} 


est borné dans Es (Qo,p). D’après la compacité de l'injection de Le (Qo,p) dans L?2(Qo,p), 
il existe une sous-suite {/(ux)}®, convergente dans L?(Qop). De (4.22), joint à l’inéga- 
lité 0 < v1 < J(v)(x) < v2 < oc, la suite EUR est convergente dans L2(0, b), ce qui 


démontre que l'injection de H1(0, b) dans L2(0, b) est compacte. 


Rappelons aussi que, de la théorie classique bien connue et du calcul explicite ré- 


sulte immédiatement que les relations correspondantes à (4.22)-(4.23) et à l’affirma- 
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tion du lemme 4.3.1 pour les espaces L?(a, b) et H'(a, b). Plus précisément, on a 


b 
[ JU(p) dx = 47 h r*|g(r)ar, (4.24) 
a,b (4) 


2 bd : 
Î JIVJ(@)(x)| dx=ar | r pu) dr, (4.25) 
Bb a dr 


et l'injection de H'(a, b) dans L?(a, b) est compacte. 
Comme l'injection de H}(0, b) dans L2(0, b) est compacte, il existe une base ortho- 
normale {e,}®, de L2(0,b) qui est un système orthogonal dans H}(0, b). Analogue- 


ment il existe une base orthonormale {dx}, de L? (a, b) qui est un système orthogo- 


nal dans A*(a, b). Les propriétés de ces bases orthonormales jouent le rôle crucial. Or, 


2 


pour démontrer les propriétés de la base orthonormale {e-}®, de L2(0,b) que nous 


=] 


allons utiliser, nous avons besoin d’une élaboration non indifférente. Donc, en ren- 


voyant à l'étude de {e,}®, à la section suivante, nous allons ici spécifier la base ortho- 


CO 


normale {d4}®., de L? (a, b) qui est un système orthogonal dans H'(a, b). 


Lemme 4.3.2. On pose 


V2 sin(ÆE(r-a)) 


dy(r) = 
Le Vb-a r 


pour a=<r=<b. (4.26) 


(A) la famille {d}®. est une base orthonormale de L?(a, b) et un système orthogonal 
k=1 


de H\(a, b) ; 


(B) on a 
ère 
re 
’ (b-— a) 


Démonstration. Comme l'injection de H'(a, b) dans L?(a, b) est compacte, il existe 
une base orthonormale {d}}®, de L? (a, b) (qui est un système orthogonal complet de 


H(a, b)) et une suite de nombres strictement positifs {Ak}. telles que 


(dk, P) EH (a,b) = Ak(di D) (ab VE H'(a, b), 
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ce qui, en vertu de (4.18), (4.20), implique que 


Prd sd b . 
- (St dm)}etnar = ax [ r'dx(r)g(r)dr Ve H'(a,b), 
où 


d 24 2 
—— (7 — dy(r)) = Azrdr(r). 

2 Fe k(r)) kr” dg(r) 

Comme on le constate sans difficulté, cette équation, jointe aux conditions d£(a) = 
_  _ : A) 
dy(b) = 0 et Idxl52ça,b) = 1, n’admet que les solutions dy(r) = = avec 
Ag=,k=1,2,.... 
D'après la définition (4.26) de d}(r) etla définition (4.20) du produit scalaire (:,)5244m 


on a 


RE kr mT 
(dk: dm)2(ub) = —| sin(o tr _ a)}sin(s = = a)\dr =Ô En 


ce qui signifie que {d}}® , est un système orthonormal de L?(a, b). 


D'autre part, comme dy(a) = dy(b) = 0 et 


_d CR) jet a) 


dr / dr r b-a r 


on a 


2 f' (4 sn (0 2))|| d re) 


(dk, dm) fi (ab) = Da r de £ dr £ 


_ 2 1H d 4 sn(pau- à) (Eat a), 


7 b-a Le r r 


(7 


ET) Oum 


ce qui signifie que {d4}® , est un système orthogonal dans H'(a, b) et 


k2 2 
CAPES 
(a,b) (b— a)? 
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On remarque aussi 


(C) on a 


Vrefa,b], VkeN, 


Idx(r)l < V2 
avb 


ce qui découle immédiatement de la définition (4.26). 


4.4 Une variante de série de Fourier sur le segment radial 


Pour construire la base orthonormale {e-}®, de L2 (0, b) qui est un système ortho- 
gonal dans À}(0, b), nous allons construire d’abord, pour chaque À > 0, les fonctions 
y(r) qui satisfont à l'équation 

d d 
(x yu) = —Av(r)r{y(r) (4.27) 


dr 


et la condition 


y(b) =0, y/() 0. (4.28) 


Par la commodité du calcul, nous posons 
s=b-r, V(s) = v(b—s), K(S)=K(b—Ss), 
et considérons l'équation 
_ [b-92r(9 Sa 9)=-AV (9-9 220;5) (4.29) 
et les conditions initiales 


d 
z(4;0) = 0, — Z(À; 5) =]; (4.30) 
ds s=0 
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Lemme 4.4.1. Le probème (4.29)-(4.30) admet une solution z(;5) et une seule sur 


l'intervalle [0, bf et elle a la forme 


Te no es b 
= D ———— pour 0<s<b-a, 
200 JAN nr (4.31) 
sin(8)+ FL (s-b+a)) 
EU ——_— pour b-a<s<b, 


où B(À) et y(À) sont définis selon les cas comme suit : 


b-a [Av 
TT K2 


VoK2 AV) K2 —K] 
À) = t t b ; 4.32 
BU) = arcco | a (1/ = (b— a)) + R) (4.32) 


sin(/#2(4-a)) 


b 
Av sin (B(À)) 


i) si2m < <2m+lavecunmeN, 


VU ; (4.33) 


ii) si2m+l< 22/2 <2(m+1) avec un mEeN, 


BQ) comme dans le cas i), 


&  Sin( 22 (b— a)) 


dd 
100) Av sin (B(À)) 


1.1 


iii) si ba 72 = 2m avec un m € N\{0}, 


(4.34) 


: . b— À 
iv) si ba > =2m+lavecunmen, 


BA =0,  yd=- 


Démonstration. Par des calculs directs, on constate que la fonction Z({; s) définie sur 
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[0, b-— a] satisfait aux conditions (4.30) et vérifie l’équation (4.29) dans l'intervalle ]0, b— 
al. 


On rappelle que, quel que soit (B, y) € R°, la fonction 


sin(B+ M(s—b+a)) 


or b-s 


satisfait à l'équation (4.29) sur l'intervalle ]b—a, b[. Or, la continuité de z({; s) et celle de 


K(S) Lz(; s) au point s = b— a impliquent que B et y doivent satisfaire aux conditions 
sin f K2 sin ( 2 (b — a) 
y =] ——b 2 , (4.35) 
a AVo a 


Av] cos sin 
K1Y P + K1Y Ê = 
K] a a 


cos | bre) G ml #2(b- a) 
—— —" " _— —————— K É eee 


a AV» a? 


= K)b (4.36) 


Pour déterminer B = B(A), y = y(À) pour chaque 1 > 0, il faut distinguer les quatre 


cas : i) sin(/22(b- a)) > 0, ii) sin(/22(b- a)) < 0, iii) sin(/22(b- a))=0et 
cos | 22 (b- a)) = 1, iv) sin( 22 (b- a)) = 0 et cos( 22 (b- a)) = —1. Ces quatre 
cas, comme on le voit facilement, ne sont que les quatre cas cités dans l'énoncé du 


lemme. Dans les cas i) et ii), en divisant les deux membres de (4.36) par 


- À 
; Av: sin 6 2 K1K2V] sin (= a)) 
K] a V2 a L 


__ fv2k2 Av2 : Ko —K] 
cot(B(A)) = a | n (b es (4.37) 


d’où l'expression (4.32)-(4.33) et celle du cas ii). Dans les cas iii) et iv), en substituant 


on obtient 


les valeurs de sin ( . (b-— a)) et de cos . (b— a)) dans (4.35)-(4.36), on obtient les 
expressions indiquées dans l'énoncé du lemme. On constate sans difficulté que dans 


tous les cas i), ii), iii), iv), la fonction z(J; 5) avec B = B(À), y = y() ainsi déterminés 
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vérifie l'équation (4.29) et les conditions (4.30). 


L'unicité de la solution z(J; 5) résulte de la théorie classique des équations différen- 


tielles ordinaires. 


Pour construire la base orthonormale {ek}®, de L2(0,b) qui est un système or- 
thogonal dans H}(0, b), certes, nous utilisons les fonctions z(Â; s) construites dans le 
lemme 4.4.1, mais il nous faut en sélectionner celles qui peuvent être un des éléments 


de cette base orthonormale. Pour cela, nous rappelons d’abord des propriétés de z(J; 5). 


Lemme 4.4.2. La fonction z(À;s) donnée dans le lemme 4.4.1 appartient à L2(0,b)n 


H (0, b) si et seulement si 


poy+ 1] Pa- mr avecunmenN. (4.38) 
1 


En outre, on a 


ion |[Z(A; s)| < co, (4.39) 
S— 
d 
—Z({; s)|5=p = 0, (4.40) 
ds 


si et seulement si (4.38) est vérifiée. 


Démonstration. Si (4.38) est vérifiée, alors, comme sin(mx -— r') = +sin(r') (le signe 


étant + ou — selon la parité de m), en posant 


on a 


ZA; 5) = +y(À) pourO<r= a. (4.41) 


On a donc 


: Avi 
[sin(/2r)| 1 
im [20:91 = lim |A 2 y] LT < oo 
s—b r—0* r K] 
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En outre, on a 
Avi 
K1 


sin( r) 


d 
51e = +7) 
as Sls=r +y( 7 


Cette dernière relation implique que X(s) _ Z(; 5) est bornée, ce qui implique que z(;:) € 
H}(0, b). 
D'autre part, si + (BG) + 4/ en) gN, alors il existe un € €]0, xl tel que 


Û Avi 
SinN(E + 1/ en 
r 


Z(4; 5) = +y() 


pour 0<r<a(b-a<s< pb). 


Dans ce cas on a 


Avi Û Av: 
d Av. CoS(e+1/ 7) sin(e+4/r) 
En. +rw| VI K] : K] | 
ds | K] r 


r2 


et donc il existe un Ô > 0 tel que 
d 1 i 
zGb- | 2 MOIS pou tout r e]0,ô[. 


On a donc 


d d 
.NI2 _ à DH Lar— 
IC; on = / KO ZA bn dr = +00. 


Maintenant nous allons chercher à caractériser les valeurs de À pour lesquelles 
z(;-) appartiennent à L2(0,b) n H}(0,b) et à montrer les propriétés utiles des fonc- 
tions y(;r) = z(1; b-r). Pour cela, rappelons que pour 0 < r < a la fonction y(A;r) = 


z(G; b-r) a l'expression (4.41) 


y@r}= 204; b-7r) =+y(Q) 


Nous commençons par le lemme suivant, qui s'écoule immédiatement d’un calcul élé- 


mentaire. 
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Lemme 4.4.3. Soit € > 0. Posons 
l—E 
Polx) = cot x — En xEeR:\(N), (4.42) 


GN=1{xeR|x = kr, k e N}). Alors @o(x) est strictement décroissante dans chaque in- 


tervalle ]kn,(k+1)r[, keN, etona 


lim @o(x) = +co, lim  @o(x) = -œ pout tout ke N. (4.43) 
x—(kr)* x—((Kk+1)7)— 


Démonstration. On a 


é (x) = -1 coba = 
270 : x2 sin?x x? 


Or, comme 


sin x < X Vx>0O, 


on a 


1 
> Fo D VxeR.,\(rN). 
sin x ZX 


poto an 


Donc @o(x) est strictement décroissante dans chaque intervalle ]kr,(k+1)7[, ke N. 


Les relations (4.43) résultent immédiatement du calcul direct. 


Considérons maintenant les deux suites {p}}#_, et {p;s}%_, définies par 
VK VK 
ha, Dép ee, (4.44) 
ay/Vi (b— a),/v2 
Posons 
E = {ppp E° = {phoh (4.45) 


en tant qu'ensembles. En ordonnant tous les nombres de E* U E** selon la grandeur 


naturelle, en comptant deux fois dans le cas où p}, = p}; et en ajoutant po = 0, on 
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CO 


définit la suite {pk}®,. On aura 


{pli = {OIUE"UE"", (4.46) 
Prk-1 <Pk<Pkr  SiPk=PpEeE"\E", 
Pk-1 <Pk<Pkr1 Sipk=Pyr € E°°\E", 


Pr-1 < Pk= Péri < Pk+2  SiPk= Péri = Ph = Pr EE NET. 


On définit la fonction ®(q) par 


V1 
®(q) = cot(a/ 219) + V 
K] 


K2V2 V2 K2—K] 

cot((b-— a); )+ (4.47) 
K1V1 ayK1v1q 
pour 


Pk-1 <q <PKk (avec prK-1 < pk). 


Lemme 4.4.4. Pout tout k € N\{0} tel que pr-1 < pr, il existe une valeur q = q+ et une 


seule appartenant à l'intervalle ]px-1, pl et vérifiant l'équation 


d(q) = 0. 


Démonstration. En posant x = 4;/ LL on réécrit la fonction D(gq) comme 
®(q) = polx) +cot(Ax), 


où po(x) est la fonction définie dans (4.42) avec € = . et 


4 (b— a)/vixki 


aw/ViK2 
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Posons 
- u A 
= KT, His KE 
On rappelle que 
lim cot(Ax) = +oo 
x—0* 
lim cot(Ax) = +00, lim cot(Ax)=-00, pour k”=1,2,3:- 
2 (int x (in) 


KX —X%X 


et que cot(Ax) est continue et strictement décroissante dans chaque intervalle ]%}», Xp, 1! 
pour tout £” € N. Ceci, joint à l'affirmation du lemme 4.4.3, implique que Q(q) est 
continue et strictement décroissante dans chaque intervelle ]p4-_:, pKl et que Q(q) vé- 


rifie les relations 


lim Q(q) = +, lim Q(q) = +0, lim_Q(q) = —-co. 
q—0* q—P£ TPE 


On en déduit l'affirmation du lemme. 


Le lemme 4.4.4 étant démontré, définissons la suite {q4}®, par les relations : 


Si Pr-1 < Ph alors px-1 << 4x < Pr Q(qy) = 0, (4.48) 


Si De-1 = Ph alors Px-1 = 4x = Pr. (4.49) 


On remarque que d’après la définition de p4, k = 0,1,2,--- (voir (4.44), (4.45), (4.46)) on 
a 


0<qg<--<qr<Gyy <<", 4x; — © pour k — co. (4.50) 


On pose 


Ak= 4%. (4.51) 
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Lemme 4.4.5. Pour les fonctions z(]; s) définies dans (4.31) on a 
lim |z(4; s)| < co, 
s— b= 


si et seulement si À = À; pour un k e N\{0}, où Ag, k = 1,2,---, sont les nombres définis 


dans (4.51). La fonction z(A y; s) possède k — 1 zéros dans l'intervalle ]0, bf. 


Démonstration. D'après le lemme 4.4.2, si lim,_.p-|Z(4; s)| < co, alors z(4; 5) = z(4; b — 


r) a l'expression 


ZA; b- 7) = +} ———— pour 0<r=<a. 
r 
Or, comme Z(J; 5) ainsi que RL z(; s) doivent être continues, on a 
lim |Z(4; 5)| <oo 
s—b7 


si et seulement s’il existe un y e R tel que 


sin( 214) R.sin(/#2(b- a) .. 
Y à — FEA RE SE , (4.52) 
d Sin Ar) NES d Sin 225) 
=ÿ= \/ Sn 4.53 
Fa P Îr=a K2 Av ds Fes Îs=b-a ( ) 


Si sin( a) =0etsin( #2(b — à)) Z 0, ou bien sin( 22 (b — a)) =0etsin( a) £ 


0, alors on voit immédiatement qu’il n'existe pas y qui vérifie l'égalité (4.52). 


sin({ M sin({/ 26 
K] K2 


D'autre part, si 
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alors (4.52) est automatiquement vérifiée, tandis que (4.53) est réduite à 


Av, cos 1 a) cos ( 22 (b- a)) 
KR ——— = Ro b———, 
K] a a 


qui nous permet de déterminer y qui satisfait aux conditions (4.52) et (4.53). 


14) £ 0 et sin(1/2*2(b— a)) £ 0, de manière analogue à la dé- 


Dans le cas sin( era . 


duction de (4.37) des égalités (4.35)-(4.36), on obtient la relation : il existeun yEeR 
vérifiant (4.52) et (4.53) si et seulement si (VA) = 0, où ®(-) est la fonction définie 
dans (4.47). 


Des considération mentionnées ci-dessus et du lemme 4.4.4 on déduit que 
lim |2(4;b-r)| <o 
r—0* 


si et seulement si À = À£ avec un & € N\{0}. En outre la définition de {pk}, et de 


{gr}, on déduit sans difficulté que y(Ax; r) = z(4; b-r) a k—1 zéros dans l'intervalle 


10, b[. 


Proposition 4.4.1. On pose 


eg(r) = nie (4.54) 
1x0) 1172ç0,m 
où 
yelr) = zQxS), (4.55) 


Z(; 5) étant la fonction définie dans (4.31). 
Alors 
(A) {ex} ®. est une base orthonormale de L2(0, b) ; 
(B) tek}, est un système orthogonal complet de H} (0, b) ; 
(EC) on a 


2 ie 
lex = #5 
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(D) eg(r) et r2x(r) Le(r) sont continues ; 
Œ) on a Æer(r)lr=0 = 0, e(b) =0; 


(F) on a 


d d 
et = Axe ext) Vrel0,aluJa,bl; 


(G) il existe une constante C\ < co telle que 


ler(a)l| < Ci VkeN\{0}, 


(H) er) possède k — 1 zéros sur ]0, bl; 


(1) on a 


VK1K2 |: 
a(a/Vik2+(b-a)/V2K)/ 


Motk-1)2<Ax<Motk+1)?,  Mo= | 
4.5 Démonstration de la proposition 4.4.1 


La démonstration de la proposition 4.4.1 s'articule en plusieurs étapes. Nous pro- 
cédons donc en démontrant d’abord des lemmes et ensuite en démontrant la propo- 


sition. 


Lemme 4.5.1. Il existe un opérateur linéaire A de L2(0, b) sur H}(0, b) qui à chaque 


u € L?(0, b) associe Au € H}(0, b) vérifiant la relation 


(Au, P) f1(0,b) = (u, P)72(0,b) VpE ete (0, b). (4.56) 
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L'opérateur À, considéré comme opérateur de L2(0, b) dans lui-même, est un opérateur 


linéaire auto-adjoint et compact et ses vecteurs propres êx correspondants aux valeurs 


propres Uk, 


AC; — er, (4.57) 


forment une base orthonormale {84}, de L2 (0, b). Les valeurs propres 1x vérifient la 
relation 


HE D DORE N 2 0 ux—0 pour k— co. (4.58) 


Démonstration. Comme on a 


b b 
Î vnréigtniar= | v(r)r di 2p0r)-Ÿ, -p(r)dr'dr 
0 0 


< cf voriene) (freres I Lot) 


avec 


on a 


Ip17200, b) — < Clplz, b) VE H(0, b). 


De la sorte le produit scalaire (u,p)5240,n Peut être considéré comme fonctionnelle 
linéaire sur @ € H}(0,b). Donc d’après le théorème de représentation de Riesz, pour 


chaque u € L?(0, b) il existe un élément U € H}(0, b) et un seul qui vérifie la relation 
{U, P) H: 0,5 = {u, P?T2(0,b) Vo E H}(0, b). (4.59) 


En posant Au = U, on définit l'opérateur A de L2(0, b) dans H}(0, b). 


Cela étant, le lemme résulte de la théorie bien connue sur l'opérateur linéaire auto- 
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adjoint et compact et de la définition des espaces L2(0, b) et H}(0, b) (voir par exemple 


le lemme 1 du $ 1 du chap. IV de [28]). 


Lemme 4.5.2. Soit {,}®, la base orthonormale de L2(0, b) définie dans le lemme 


4.5.1. Alors pour tout k on a 


df od_). . : 
per —&k) sn. dans L°(0, b). (4.60) 


Démonstration. Comme Aëk = uKe4, en posant U = A6y, d’après (4.59) on a 


(La) Letar = [vor 
Lux (a) puar= EURO ek(r)p(r)dr. 


Compte tenu que A} (0, b) est dense dans L2?(0, b), cette égalité définit la dérivée géné- 


ralisée L(x(r) r2-La(r)), qui vérifie (4.60). 


Lemme 4.5.3. Soit {&,}®, la base orthonormale de L2(0, b) définie dans le lemme 


4.5.1. Alors pour tout k, ë(r) a la forme 


> yk(r) 


A » (4.61) 
I yx(r) 1172(0,b) 


où yK(r) est la fonction définie dans (4.55). 


Démonstration. D'après le lemme 4.5.2, e4(r) doit vérifier l'équation (4.60) qui est iden- 
tique à l'équation (4.27) avec À = . et vérifier aussi les conditions aux limites (4.28). 
Donc &;(r) doit avoir la forme de la fonction y(4;r) = z(4; b-r) (définie dans (4.31)) 
avec À = rs Or, êg(r) doit appartenir à L2(0, b) n H}(0,b) et donc, d’après le lemme 


4.4.5, il faut et il suffit que À = ma soit une des valeurs de {A4}® .. Comme {uk}, est 
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ordonnée dans l’ordre décroissant tandis que {AH}? . est ordonnée dans l’ordre crois- 
sant,ona 


À = —. (4.62) 


C'est-à-dire, &-(r) a la forme (4.61). 


Maintenant nous démontrons la proposition 4.4.1. 


Démonstration. Comme les relations 


sont établies (le lemme 4.5.3), la relation (A) est démontrée. La relation (B) résulte im- 
médiatement du lemme 4.5.1. Pour établir (C) il suffit de substituer u = @ = ex dans 
(4.56) et de tenir compte de (4.57). Les relations (D) et (F) résultent immédiatement de 
(4.60). D'autre part, les relations (E) et (H) sont établies dans la définition de e-(r) et 
dans le lemme 4.42. 

Pour démontrer l’inégalié 


lek(a)l < Ci 


il nous est commode de considérer la fonction 


: At 
K2 ee V à 


AgV2 b-s 


zZ( 4; 5) = pour 0<s<b-a, 


de sorte que 


sin ( ? 
SR PR Re 
AgV2 a 


Az 


eg(a) = 


De cette dernière expression on obtient immédiatement 


b l 
lex(a)| < - i 
Vavalzd)lzon 
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Considérons d’abord le cas où 


c'est-à-dire le cas où 


Dans ce cas on a 


b-a 
sine) £ $) = 
ee 9], RE 
= #20 082 [b Ge ni Dee [= 
2 V Agv 2/AgVv2 K2 4 V Agv 
vob K2 1 zx K2 
> 1 b 1/ > 
É 2 AgV2 | S G ü à nl . 


v2b K2 
SE En ET 
On en déduit que 
F0 blGr-b@-avb E-av 


Comme l’ensemble des nombres À} tels que 


T°K2 


À LR Er ee, 
ES peas 
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est fini (et que e&(r) sont continues), on a 


sup |eg(a)|= (ele < O0 
keAo 


où 


H°K2 


Â0 = {RENE < ———"—}. 
o—{ x bar 
On en déduit que 


sup lex(a)| < max{Cl*!, ea = C1 <oo. 
ke Ao 


De l'expression de eg(r) (voir (4.54), (4.55), (4.31)) on déduit que sur l'intervalle 


10, a] lya[VA Se 


entière d’un nombre x > 0). Comme le nombre de zéros de e-(r) sur ]0, b[ est k—1,on 


[VAx = PRE x] zéros (ici [x] est la partie 


a 
ay (b-— a) 
k=| À” 4e +| Re] 
2 (b—a)ÿv2 (b— a),/v2 
er LÀ (et Vi V4 xp 


et 


k= (VHS) - RE) +1 autrement. 


Comme pour tout (a, B) ER; x*R;ona 
a+B-2<[a]+[fl<a+f, 


on a 


ju Ne (b— a), /vV2K] on | Nr ARR DAVVERT, 


VK1K2 _ _ VK1K2 
d'où 
VRE — 1) VRTRZ 


< 


RE — — SE 
T(a1/ViKk2 + (b—4a)/V2K1) ” x(ay/Vik2 +(b- a)/V2Kk1) 


= 


(k+ 1). 
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La propriété (1) est démontrée. 


4.6 Solution stationnaire 


Maintenant nous allons chercher la solution stationnaire (T(r),II(r)) de notre pro- 


blème. Plus précisément, nous allons construire les fonctions T(r) et [I(r) qui satisfont 


aux équations 


nl 


avec y définie par 


[| 
Y = Yi 


T=Q+E) [| r=a 


€] 
et les conditions aux limites 


NÉcSeTs, 
IT |;=4= Ta +] T |r=a 
Lies 


On définit en outre 6, par l'égalité 


. A 
O= Te (To - Ga + (To + Ou), 
c'est-à-dire 
: A : 
GP qu 
b+7y1A4@1 


AD HO. 


On a la proposition suivante. 


xD ET) =wyô(r-a) pour 0<r<b, 


d 
(F0) = 0 pour a<r<b 


(4.63) 


(4.64) 


(4.65) 


(4.66) 


(4.67) 


(4.68) 


(4.69) 
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Proposition 4.6.1. La solution du problème (4.63)—(4.68) est donnée par 


_ T1, (Gr, -aill, 10 .)er(a) 
Tr) = Tp+ D et, =, 
le=il b À& 


bTip — Aa + & (To + Oa)) … Ab(Iy — (Ta + (Ty + Oa))) 


Los (b=a b=ar 


(4.70) 


(4.71) 


Démonstration. Il est clair que la solution (formelle) du problème (4.64), (4.67), (4.68) 


est 
by — aa + T(a)) _ ab(y - (Ta + M T(a))) 


Hu) = (b-a b=ar 


Donc, en substituant cette relation dans (4.65), on obtient 


de. | =. b(Tl — (Ta + 1 T(a))) 
l l 


Per €] : (b— a)(a+E€1) 


Nous posons également 


O(r)= T(r)- Th. 


Alors le problème (4.63), (4.66) se réduit à 


4.5.4 _ l 
et «(r) 2 @(r) = Ya(T(a))ô(r — a), 


O(b) = 0. 
En posant 


@(r) = Y cPer(r) 
k=1 


= Ye, (T(a)). 


(4.72) 


(4.73) 


(4.74) 


(4.75) 


(4.76) 


et en intégrant les deux membres de (4.75) multipliés par e(r) sur ]0, b], en vertu de 


(B), (C), (D) et (F) de la proposition 4.4.1,ona 


2 d d 
Me = fetn(-(rktntn)}ar = vatriaerta 
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C'est-à-dire, la solution du problème (4.75)-(4.76) est donnée par 


eu = À Cet, = Yatrta) À. (4.77) 
k=1 À 
Donc la valeur de @(a) est 
CO À 
E(@ = aa) 4e  Aa= 9 EN, (4.78) 
k=1 À 
De (4.73) et de (4.78) on obtient 
AD 
(a) = —] (Tip — (Ta + 1 (T + O(Q)))). (4.79) 


(b— a)(a+€1) 


Donc, d’après la définition (4.69) de 6,, on a 


O(a) = 84 


et donc 


T(a) = Ty +4. 


En substituant cette expression dans (4.73), (4.77), (4.72), on obtient (4.70) et (4.71). 


4.7 Cas du domaine d’une dimension spatiale 


Dans la section précédente nous avons déjà montré l'existence et l’unicité de la 


solution stationnaire de ce problème, que nous notons (Tr), I,:(r)). Posons 
d=0(t,r)=T(t,r)-Tr), (4.80) 


n=ntr)=0(Gr)-ss(r) - 0 (47). (4.81) 
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Comme 6,(r20,I1) = 0 (voir (4.64)), de (4.2) découle 
On Yoôr(r/ôrn) = -&1(0:9—Yoôr(r/0,9)) pour a<r<b. 


On remarque que, en vertu de (4.1), (4.4),(4.5) et (4.80), on a 


_ K@) 
—%1 (09 — yoôr(r/0,8)) = n(— x — o]or(r 20,8) pour 0<z< a. 
Cvo 


Ces relations nous permettent de réduire le système d'équations (4.1)-(4.3) à 


Cup0r9—0,(r"k0,9) = qô(r-a) dans ]0,bI, (4.82) 
0sn—Yyoô-r(r”ô;n)= f dans Ja,bl, (4.83) 
q = q{) ee (4.84) 
€] 
= a(S- roJôr(r?0,8). (4.85) 
Cve 


Les conditions aux limites sont 
O(E, b) =0, n(t, a) = n(t,b) =0 Vt2>0, (4.86) 
tandis que les conditions initiales sont 
9(0,r) = To(r)-T tr) = 8o(r) pour0<r<b, (4.87) 


n(0,r) =Ilo(r)—Ilss(r)— &100(r) =notr)  poura<r<b. (4.88) 


Pour les fonctions 4 définies sur ]0, b], en utilisant les fonctions e4 définies dans 


(4.54)-(4.55), on définit les coefficients de Fourier &(k) et la norme ||u|| He, Par 


A = (er ob ll, = (2 ARE IE (4.89) 
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où ne (0, b) est défini comme dans (4.15). D'autre part, pour les fonctions u définies 


sur {a, b], en utilisant les fonctions d4(r) définies dans (4.26) on définit 


AE) = (u,dé) ça ul, PLU (4.90) 


Proposition 4.7.1. Soit 1 < o < à. Supposons que do € H(, “ etno € He él Alors, quel 
que soit t > O, il existe une solution (9,1) et une seule du problème (4.82)-(4.88) dans 
la classe 


1 2(n +. 
ÿe L°(0, 6 Ho a) NL (0,5 HG) 


DEL (0H DM, 4) (4.91) 


Démonstration. En utilisantles coefficients de Fourier définis dans (4.89)-(4.90) (même 
si on utilise le même symbole #(k) pour les deux sens différents, il n’est pas difficile de 


les distinguer du contexte), on peut transformer les équations (4.82)-(4.83) en 
d = 2 
at k) + ALÜ(E, k) = q(Per(a), (4.92) 


d : 
CD + KA = FD). (4.93) 


a 


En multipliant les deux membres de (4.92) par 19-18 (4, k),ona 
1 = = ne 
AA LE, k)) + ALTO (1, © = AT 9 (4, k)q(der(a). 


Comme 
lex(a)| < Ci VkeN\{0} 
(voir (G) de la proposition 4.4.1), on a 


2 > CF 
ALTO, )q(Der(a) < ALIÈE, HI + SAR NO. 
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Or, comme o —2 < — d en vertu de (1) de la proposition 4.4.1 on a 


OO 
à À = Ci <oo. 
k=1 


On en déduit que 


d 
alle + le, < CC1lqO. 


Quant à |g(#l*,ona 


2 2 
IS. 2 
POIL Titre = À nt Haen| = 
l 1 = 


OO 
< CL ET MG DE = Clin 
k=1 Ho 


2 1 


Æ <oo. 
4 E Pb- à À Fi L 


En adjoignant les deux inégalités, on obtient 


d 3 2 2 2 

—_— — : 4.94 
ai El no + Me CiCaCalni s-J (4.94) 
(a,b) 


D'autre part, si on multiplie les deux membres de (4.93) par k27-2f(#, k) et on fait 


la somme par rapport à k, on obtient 


d 2 T 2 2 
— 901 + Vo lnlge  < Cf (EE 0-2 » 
dt % H 5 LL (b -— a)? 1 Hop f Hp 


(4.95) 
où 
1 . 2 
ee (b . | 
Yo 7 
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Si on définit 
ÿ=e "9, G=e "in q=e "ta, f=e"!f (4.96) 
avec un & > 0 (à choisir dans la suite), de (4.82)-(4.83) on obtient 


Coo019 — 0, (r?k0, Ÿ) + @CypŸ = qô(r) dans ]0, b[, (4.97) 


05 yoôr(r/0,)+w=f dans [a,bl. (4.98) 


En procédant d’une manière tout analogue à l'obtention de (4.94)-(4.95), de (4.97)- 


(4.98) on obtient 


TIGE 1 La lle, +201 1 = CCaCqlAle ne: (4.99) 
HG, b) 
d 
di Mae TER re +20 re de < Cl flr- 2° (4.100) 
Soit 9 une fonction appartenant à la classe 
eo eh nt 4) (4.101) 


Alors, en substituant ® à la place de ® dans la définition (4.85) de f, on définit f, qui 
appartiendra à L?(0, f; H (al En substituant f à la place de f dans (4.83), on peut ré- 


soudre l'équation (4.83) (avec f) et on obtient 
EL, 6 H7 NL (0,5 4H) 


(voir (4.95)). Ensuite, en définissant gq par (4.84), on résout l’équation (4.82) et on ob- 
tient 0, qui, en vertu de (4.94), apprtient à la classe Y. C'est-à-dire, on a définit un 
opérateur G1 qui, à 9 € Ys, associe la solution 9 € Y3 de l'équation (4.82). 


Si on définit f = e “! f et si on rappelle la définition de f, on voit aisément qu’il 
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existe une constante C; telle que 


Ie & Cr : (4.102) 


Considérons maintenant deux fonctions À et % appartenant à la classe Yo. En 


définissant #1, 0, :, Do, T1, 2 de la manière indiquée ci-dessus, on pose 
Ô =, À = 1 — 2, O = 91 — 0. 


Comme les équations (4.97)-(4.98) sont linéaires, de la même manière que (4.99)- 


(4.100) et (4.102) on obtient 


TION + Ole, + 2016 + à < CIC A, sis (4.103) 
Le 
D 
MAÉ 1 0 al Axe, +201 He à < CCE, ; (4.104) 
Comme 4271 < 642 + 42072 et donc ||H|? FE SU AlGe , 9 + LIAIE Ho-1 POUr 


(a,b) 
tout Ô > 0, en multipliant les deux membres de (4.103) par une constante À à ae 


ner et en faisant la somme de ces deux inégalités, on obtient 


d 
AIG ge + IA 


+ AIO + 
al I le 


ue 1) 


2 
JT _ 
og lue, + 2A@ION 0 + 2@ 


(b— a)? HD 
2 12 À > 2 2 
< SAC Ca Cal A + 3 C1 Ca Cal As à + CCF . (4.105) 
2 
Si on choisit À = 2C2Cf, ô = Aa da = À C?C1C l'inégalité (4.105) se 
réduit à 
d A 12 T 2 
PAU us l'A ge 1)+AI le, V0 a 2(b- a} IL Ale, 


Â 


on SO, 


+26 so. 1 + &| AI (4.106) 
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L'inégalité (4.106) implique que l’approximation successive Ü (ny n=1,2;,..., pour 
le problème (4.97)-(4.98) (avec les conditions initiales (4.87)-(4.88) pour tout n avec 


do € H ce et 70 € H É cn) converge dans l’espace de Banach Ÿ) muni de la norme 


c. 7 = - 1/2 
— 2 2 2 
ID, = (El one + oz) + 2010 ze) 
Comme les équations sont linéaires, la limite Ÿ de l’approximation successive Ü(» 
sera la solution du problème (4.97)-(4.98) (et, par sa construction, donnera aussi }). 
En posant 9 = e”! ë, n = e*!5, on voit aisément que (#, n) sera la solution du problème 


(4.82)-(4.88). Leur appartenance aux classes (4.91) résulte des estimations établies ci- 


dessus. 


La solution (0,7) du problème (4.82)-(41.88) étant trouvée, par les égalités (4.80)- 
(4.81) on peut construire les fonctions T(f,r) et II(f,r) qui satisfont aux équations 


(4.1)—(4.3) et aux conditions (4.6)-(4.8). 
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Ss 


AANS la présente thèse, nous avons démontré l'existence et l’unicité de la solu- 


L 


tion du système de deux équations de type parabolique linéaires, qui prend 
son origine de la modélisation du phénomène de l’évaporation de l’eau à partir de la 
surface d’eau liquide, et ce dans le cas d’un domaine d’une dimension spatiale, dans le 
cas d’un domaine de dimension 3 délimité par deux plans horizontaux ainsi que dans 
un domaine sphérique avec la symétrie sphérique, même si dans le cas du domaine de 
dimension 3 le résultat a été obtenu sous une hypothèse assez restrictive. 

Compte tenu de la difficulté intrinsèque de la modélisation mathématique du phé- 
nomène de l’évaporation (et de son effet par la chaleur latente), difficulté essentielle- 
ment due au fait que l’évaporation se produit sur la surface du liquide — variété de di- 
mension 2 -, tandis que l’on est intéressé au comportement des grandeurs physiques 
(la température et la densité de la vapeur) dans une région de dimension 3, nous pen- 
sons que le résultat obtenu est un bon indicateur du comportement de la température 
et de la densité de la vapeur dans un domaine où se produit l’évaporation, bon indica- 
teur dans la possibilité de l'étude dans le cadre de l'Analyse mathématique. 

Pour obtenir ces résultats, nous avons développé une méthode basée sur des va- 
riantes de série de Fourier. En concret, il s’agit d’une variante de série de Fourier sur 
un intervalle ] — b,a[ composé par deux sous-intervalles ] — b,0[ et ]0, al où les co- 


efficients de l’opérateur différentiel sont des constantes différentes sur chaque sous- 
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intervalle, et aussi d’une variante de série de Fourier sur le rayon d’une sphère avec 
un comportement analogue des coefficients de l’opérateur différentiel. Nous pensons 
que cette méthode pourra être ultérieurement développée pour être appliquée à divers 
problèmes d'équations du type parabolique et aussi d'équations du type elliptique. 

Certainement la question qui reste pour le moment sans être résolue, celle sur la 
possibilité de nous libérer de l'hypothèse restrictive dans le cas du domaine de dimen- 
sion 3, devra être investiguée. Il se peut que cette question puisse être résolue par une 
élaboration ultérieure de calcul, mais il se peut aussi qu’à la résolution de cette ques- 
tion puisse être servir l’approfondissement de l’utilisation des espaces de Sobolev ani- 
sotropes. Dans ce dernier cas, on aura besoin d’une étude approfondie sur les espaces 
de Sobolev anisotropes. 

Pour conclure nous citons deux thèmes de recherche, qui ne sont pas dans la portée 
immédiate du travail présenté dans la présente thèse, mais auxquels notre intérêt peut 
être dirigé. Le premier est une modélisation mathématique plus cohérente du phéno- 
mène d’évaporation de l’eau et de son effet dans le milieu environnant. Le deuxième 
est l'étude du système d'équations décrivant non seulement l’évaporation de l’eau et 
son effet thermique, mais aussi le mouvement de l’eau et de l'air. Ce sont des thèmes 
de grand intérêt scientifique et on espère que les études effectuées dans la présente 


thèse pourront donner des contributions. 
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